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COMPUTING INTEGRAL BASES VIA LOCALIZATION AND HENSEL 

LIFTING 


JANKO BOHM, WOLFRAM DECKER, SANTIAGO LAPLAGNE, AND GERHARD PFISTER 

Abstract. We present a new algorithm for computing integral bases in algebraic function 
fields, or equivalently for constructing the normalization of a plane curve. Our basic strategy 
makes use of localization and, then, completion at each singularity of the curve. In this way, 
we are reduced to finding integral bases at the branches of the singularities. To solve the latter 
task, we work with suitably truncated Puiseux expansions. In contrast to van Hoeij’s algorithm 
[21], which also relies on Puiseux expansions (but pursues a different strategy), we use Hensel’s 
lemma as a key ingredient. This allows us at some steps of the algorithm to compute factors 
corresponding to complete sets of conjugate Puiseux expansions, without actually computing 
the individual expansions. In this way, we make substantially less use of the Newton-Puiseux 
algorithm. In addition, our algorithm is inherently parallel. As a result, it outperforms in 
most cases any other algorithm known to us by far. Typical applications are the computation 
of adjoint ideals [4] and, based on this, the computation of Riemann-Roch spaces and the 
parametrization of rational curves. 


1. Introduction 

Let A be a reduced Noetherian ring, and let Q(A) be its total ring of fractions. The normal¬ 
ization of A is the integral closure of A in Q(A). We denote the normalization by A and call A 
normal if A = A. Recall that if A is a reduced affine (that is, finitely generated) algebra over a 
field, then A is a finite A-module by Emmy Noether’s finiteness theorem (see [13], [20]). 

In this paper, we are interested in the case where A is the coordinate ring of an algebraic curve 
defined over a field K of characteristic zero. More precisely, let / G K[X,Y] be an irreducible 
polynomial in two variables, let C C A^(iG) be the affine plane curve defined by /, and let 

A = K[C]=K[X,Y]/{f{X,Y)) 

be the coordinate ring of C. We write x and y for the residue classes of X and Y modulo 
/, respectively. Throughout the paper, we suppose that / is monic in Y (due to Noether 
normalization, this can always be achieved by a linear change of coordinates). Then the function 
field of C is of type 

K{C) = Q(A) = K{x)[y] = K{X)[Y]/{fiX,Y)), 

a: is a separating transcendence basis of K{C) over K, and y is integral over A[a:], with integrality 
equation f(x, y) = 0. In particular, A is integral over iL[x], which implies that A coincides with 
the integral closure K[x] of K[x] in K{C). We may, hence, represent A either by generators 
over A or by generators over A[x]. For the latter, note that A = K[x] is a free iL[x]-module of 
rank 

n ■■= degyif) = [K{C) : K{x)]. 

Indeed, this follows by applying [19, Theorem 3.3.4] to the PID 

K[x] C K{x) C K{C) = K{x)[y]. 

Definition 1.1. If i? is any ring and R is a reduced Noetherian ring such that R is a finite free 
R-module, then an integral basis for B over R is a set bi,... ,br of free generators for B over R: 

B = Rbi © • • • © Rby 
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Remark 1.2. In the context outlined above, there is always an integral basis for A over K\x\ 
of type 

pi{x,y) pn-i{x,y) 

’ d{x) ’■■■’ d{x) ’ 

with d G K[x\, and with elements pi G K[x][y] of degree z in y. Such a basis is obtained from 
any given set 1 = cq, ... , Cm-i of -module generators for A by unimodular row operations 
over the PID K[x]: For each i, write c* = with coefficients Cij G K{x). Then 

take d to be the least common denominator of the c^-, transform the matrix {d ■ Cij) into row 
echolon form {pij), and set Pn-i = ^ for i = 1,..., n — 1. 

Remark 1.3. General normalization algorithms are presented in [16], [5]. They are designed 
to return an ideal U <Z A together with an element d £ A such that A = C Q(^). Here, 
as we will recall in Section 2, any non-zero element of the Jacobian ideal M of A = K\x^y] 
can be taken to be d. In particular, we can choose d to be a generator of the elimination ideal 
Mr\K[x\. The roots of d in the algebraic closure K oi K are then precisely the x-coordinates of 
the singularities of the curve defined by / in If uq = d(x), ui,... ,Ur generate the ideal 

U, the y'^Uj{x,y)/d{x), 0<z<n — 1, 0<j<r, generate A over K[x]. An integral basis is then 
obtained by operations as described in Remark 1.2. 

Remark 1.4. In practical terms, uq, ■ ■ ■ ,Ur are given as polynomials in K[X,Y] of T-degree 
at most n — 1. If these polynomials, together with /, form a Groebner basis with respect to 
the lexicographical ordering, taking Y > X, then already the elements y^Uj{x,y)/d{x), 0 <i < 
n — 1 — deg{uj), 0 < j < r, generate A over K[x]. 

Example 1.5. Gonsider the standard cusp: Let 

A = K[x,y]=K[X,Y]/{Y^-X^). 

As a module over A, we may represent A as 

A = A- — + A-l = - {y^, x)^ 

(see [16, Example 2.5]). Gonsidering A over iL[x], we get 

_ /jy2 y.2 

A = K[x\ ■ -h K[x] ■ y ■ -h K[x\ • --h K\x\ ■ 1 J- K[x\ ■ y + K[x\ ■ y^. 

XXX 

Since y^ = and K[x\ • y^ C K[x\ ■ y"^jx, we have 

_ y^ 

A = K[x\ ■ — © K[x\ ■ 1 © K[x\ ■ y. 

Hence, l,y,y‘^lx is an integral basis as in Remark 1.2. 

The algorithms in [16], [5] work for any reduced affine algebra A over a perfect field. They 
rely on the Grauert and Remmert normalization criterion which can be applied in a global or 
local setting (see [15], [17, Prop. 3.6.5], [5, Prop. 3.3]): Whereas the algorithm in [16] is of 
global nature, the idea in [5] is to consider a finite stratihcation of the singular locus Sing(A), 
apply a local version of the normalization algorithm at each stratum, and hnd A by putting 
the resulting local contributions together. If Sing(A) is finite, we may stratify it by considering 
each P G Sing(A) separately. This applies, in particular, to the case where A = K[C] is the 
coordinate ring of a curve C as outlined above. As a consequence, computing an integral basis 
for A over K[x\ is then equivalent to computing a local contribution to A at each P. 

In this paper, we present a new method for computing the local contributions which is custom- 
made for the case A = K[C]. We proceed along the following lines. To fix our ideas, in Sections 
2 and 3, we brieffy recall the Grauert and Remmert type algorithms. Furthermore, we discuss an 
efficient criterion for detecting whether a given point is the only singularity of the curve under 
consideration. In Section 4, we review the theory of Puiseux expansions and its connection to 
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integrality. In Section 5, taking an analytic point of view, we show how to obtain an integral 
basis at a given singularity P from integral bases at the branches of the singularity. In Section 6, 
we explain how to construct the local contribution at P from the integral basis at the singularity. 
How to actually find the integral bases at the branches is a topic of Section 7: Working with 
approximations by suitably truncated Puiseux series, we describe a way of writing down an 
integral basis for a single branch without performing too many computations. This approach 
is inspired by van Hoeij’s paper [21], but pursues a different strategy, with Hensel lifting as a 
crucial new ingredient. Moreover, we modify the theoretical results of Section 5 in order to 
achieve a better performance. 

We have implemented our algorithm in the open source computer algebra system Singular 
[11]. In Section 8, we compare the performance of the algorithm with that of the local to global 
approach from [5]. We also give timings for the implementation of van Hoeij’s algorithm in 
Maple and for the variant of the Round 2 algorithm implemented in Magma. 

2. The Global Normalization Algorithm 

In this section, we review the global version of the normalization algorithm. To begin with, we 
fix our notation and give some general facts on normalization. For this, A may be any reduced 
Noetherian ring. We write 


Spec(A) = {P C A \ P prime ideal} 

for the spectrum of A. The vanishing locus of an ideal J of A in Spec(A) is the set V{J) = {P G 
Spec(A) I P D Jj. We denote by 

N{A) = {P G Spec(A) I Ap is not normal} 

the non-normal locus of A, and by 

Sing(A) = {P G Spec(A) | Ap is not regular} 

the singular locus of A. Then N{A) C Sing(A), with equality holding if A is the coordinate ring 
of a curve (see [9, Theorem 4.4.9]). 

Definition 2.1. The conductor of A is 

Ca = Ann/i(A/A) = {a G A | aA C A}. 

Note that Ca is the largest ideal of A which is also an ideal of A. 

To emphasize the role of the conductor, we note: 

Lemma 2.2. Let A be a reduced Noetherian ring. Then N{A) C V{Ca)- Furthermore, A is a 
finite A-module iff Ca contains a non-zerodivisor of A. In this case, N{A) = V{Ca)- 

Note, however, that the conductor can only be computed a posteriori when A is already 
known. 

Definition 2.3. Let A be a reduced Noetherian ring. A test ideal for A is a radical ideal 
J C A such that V{Ca) C V{J). A test pair for A consists of a test ideal J together with a 
non-zerodivisor g € J of A. 

Test pairs appear in the Grauert and Remmert normality criterion which is fundamental to 
algorithmic normalization (see [15], [17, Prop. 3.6.5]). The normalization algorithm of de Jong 
(see [8], [10]) and its improvement, the algorithm of Greuel et al. [16], are based on this criterion. 
Both algorithms apply to any reduced affine algebra A = K[Xi, ..., Xn\fl over a perfect field K. 
By means of primary decomposition, we may reduce to the case where A is equidimensional. In 
this case, since we work over a perfect field, the Jacobian ideal^ M of A is non-zero and contained 

^The Jacobian ideal M of A = K[Xi ,..., X„]// is generated by the images of the c x c minors of the Jacobian 
matrix ( J(|g), where c is the codimension, and fi,..., p are generators for I. By the Jacobian criterion, V (M) = 
Sing(A) (see [13, Theorem 16.19]). 



4 


J. BOHM, W. DECKER, S. LAPLAGNE, AND G. PFISTER 


in the conductor Ca, so that we may choose the radical J = y/M together with any non-zero 
divisor in J as a test pair (see [16, Lemma 4.1]). The idea of hnding A is then to successively 
enlarge A by finite ring extensions Ai+i = Hom^^(Ji, Jj) = -Ai Ji) C A C Q(^), with 

Aq = A and Ji = y/JAi, until the normality criterion of Grauert and Remmert allows us to 
stop. As already pointed out in Remark 1.3, the algorithm of Greuel et al. is designed so that 
it returns an ideal U <Z A together with an element d ^ A such that A = C Q(A). 

Remark 2.4. If M is non-zero and contained in Ca, then any non-zero divisor in M is valid as 
a denominator: If 0 7 ^ c G M, and A = as above, then c ■ =: t/' is an ideal of A, and 

lU=lU'. 

Example 2.5. Let A be the coordinate ring of the curve C with defining polynomial f(X,Y) = 
X5 _ g Q^x, Y]. Then 

J := (x,y(y- 1))^ 

is the radical of the Jacobian ideal, so we can take (J, x) as a test pair. In its first step, the 
normalization algorithm yields 

-4i = -Ui = ^ {x,y{y - if)^. 

X X ^ 

In the next steps, we get 

A 2 = ^U 2 = ^ (x^xy(y - l),y(y - if)^ 

and 

^3 = ^t ^3 = ^ {x^,x‘^y{y - 1 ), xy{y - if, y‘^{y - 1 )^)^ . 

In the final step, we find that A 3 is normal and, hence, equal to A. 

3. Normalization of Curves via Localization 

In this section, we discuss the local to global variant of the normalization algorithm proposed 
by Rohm et al. [5]. To simplify our presentation, we focus on the case of a reduced Noetherian 
ring with a finite singular locus (which includes our case of interest here). Our starting point 
is Proposition 3.1 below which is also fundamental to our new algorithm. In formulating the 
proposition, if P G Spec(A) and A C A' C A is an intermediate ring, we write A'p for the 
localization of A' at A \ P C A'. 

Proposition 3.1. Let A he a reduced Noetherian ring with a finite singular locus Sing(A) = 
{Pi,...,Ps}. For i = l,...,s, let an intermediate ring A C A^®) C A be given such that 
Ay = Ap.. Then 

^ A(®) = A. 

i=l 

Proof. A more general result is proved in [5, Proposition 3.2]. □ 

Definition 3.2. We call any ring A^®) as in the proposition a local contribution to A at Pi. If 
in addition Ap^ = Ap^ for j i, we speak of a minimal local contribution to A at Pi. 

Remark 3.3. Note that such a contribution is uniquely determined since, by definition, its 
localization at each P G Spec(A) is determined. 

Given a reduced affine algebra A over a perfect field K with a finite singular locus, Proposition 
3.1 allows us to split the computation of A into local tasks at the primes P, G Sing(A). One way 
of hnding the minimal local contributions A^*) is to apply the local version of the normalization 
algorithm from [5] which relies on a local variant of the Grauert and Remmert criterion. For 
each i, the basic idea is to use Pi together with a suitable element gi of the Jacobian ideal instead 
of a test pair as in Dehnition 2.3. 
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Example 3.4. As in Example 2.5, let A be the coordinate ring of the curve C with defining 
polynomial f{X, Y) = — Y‘^{Y — 1)^ G Y]. Note that C has a double point of type A 4 at 

(0,0) and a triple point of type Eg at (0,1). If we apply the strategy above, taking Pi = (x, 

P 2 = {y — 1, x) and gi = 02 = x, we get local contributions J-Ei, i = 1,2. Specifically, 

di=x^ and Pi = (x^, y(y - 1 ) 3 )^, 

d 2 = x^ and U 2 = (^x^,x‘^y'^ {y - 1) ,y‘^ {y - 1)'^'^ 

Summing up the local contributions, we get A = with d = and 

U = (^x^, y{y - l)^x, (y - 1) x^ y^ (y - 1)^^^ . 

Note that U coincides with the ideal P 3 computed in Example 2.5. 

Remark 3.5. In Example 3.4, the normalization of the local ring Ap^ is Apj = ^{x^,x‘^{y — 
l),(y — Indeed, since y^ is a unit in Apj, this follows by localizing U 2 at P 2 . Note, 

however, that (y—l)/x and (y—l)^/x3 are not integral over A. Hence, ^(x^, x^(y—1), (y —1 )^)a 
is not a local contribution to A at P 2 - 

Relying on the Jacobian criterion, we may find the primes in Sing(A) by means of primary 
decomposition. If there is precisely one such prime, this requires (possibly expensive) compu¬ 
tations which are only needed to detect this fact. In the case of a plane curve C considered 
here, supposing that one singularity P of C is already known to us, we may check whether P is 
the only singularity of C by comparing the local Tjurina number of C at P with the total Tju- 
rina number of C. Computing the total Tjurina number via Grobner bases over the rationals, 
however, can be expensive due to coefficient swell. To overcome this problem, we provide an 
efficient modular criterion. Note that though singularities at infinity do not matter for obtaining 
integral bases, the criterion takes these singularities into account. That is, it is formulated in 
the projective setting. 

Let K be any field, let F G K[X, Y, Z] be a square-free homogeneous polynomial of positive 
degree, and let T = Proj(iL[A, Y, Z]/{F)) be the projective curve defined by F. Moreover, write 

S = K[X,Y,Z]/{Fx,Fy,Fz), 

where Fx,Fy,Fz are the partial derivatives of P. Then, taking Euler’s rule into account, 
Tsi ng = Proj(5) C T is the singular locus of P. For any Q G Tging, let be the homogeneous 
localization of S' at Q. Then 

XQ(r) = dim^ S(Q) 

is the Tjurina number of T at Q. For example, if P = {X, Y), then 

tp(F) = dimi^ {K[X,Y\^x,Y)/{f, fxjv)) , 
with / = F{X, Y, 1). The Tjurina number of F in the chart A 7 ^ 0 is 

Tx^o(r) = dimx {K[X, Y]l if, fx, fv )), 
and similarly for the other coordinate charts. Finally, 

r(r) = degProj(S) = ^ tq 

Q^^sing 

is the total Tjurina number of F. 

Proposition 3.6. Let F G Q[A, Y, Z\ be a square-free homogeneous polynomial of positive degree 
with integer coefficients. Let q be a prime number such that the reduction Fq of F modulo q is 
non-zero. Consider the curves F = Proj(Q[A, T, E]/(P)) and Fg = Vio]{¥q[X,Y, Z]/{Fq)), and 
let P = {X, Y). Suppose that 

xp{Tq) = rp(F) > 0 and TxjLo{Tq) = TY^o{Tq) = 0. 


Then Fsing = {P}. 
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Proof. By [2, Theorem 5.3], considering the Hilbert functions of 

S = Q[X,Y,Z]/{Fx,Fy,Fz) and 
Sg = F,[X, y, Z]/ {{Fx)„ {Fy\, (Fz),) , 

we have 

HFs(t) < HFs^(t), for all t. 

Since the Tjurina numbers are the leading coefficients of the respective Hilbert polynomials, this 
implies that 

^(r)<T(F,). 

On the other hand, if Txj^ofPq) = Ty^o(rg) = 0, then (Fq)sing = {T”}, so that 

r(r,) = rp(F,) = rp(F) < t(F). 

Combining both inequalities yields 

rp(r) = t(F) 

and, thus, Fsmg = {P}. □ 

Remark 3.7. The invariants in the criterion can be obtained efficiently by a standard basis 
computation over Q with respect to a local ordering and by standard basis computations over 
Fp with respect to a global and a local ordering, respectively. 


4. PuiSEUx Series and Integrality 

We discuss some basic facts about Puiseux series and their connection to integrality. 

4.1. Puiseux Series. Let K G L he a field extension, with L algebraically closed. The field of 
Puiseux series over L is the field 

OO 

L{{X]] = U L((X‘/”)). 

m=l 

The Newton-Puiseux theorem, which is closely related to the aforementioned finiteness theorem 
of Emmy Noether, says that L{{X}} is the algebraic closure of L{(X)). In particular, L[[X^/™']] 
is the integral closure of L[[^]] in L((X^/™')). See [13, Chapter 13], [1, Lecture 12]. 

We have a canonical valuation map 

V : L{{X}} \ {0} Q, 7 hA ^(y), 

where v{'y) is the smallest exponent appearing in a term of 7 . By convention, u(0) = 00 . 
The corresponding valuation ring L{{X}}^>o consists of all Puiseux series with non-negative 
exponents only. Henceforth it will be denoted by Vx- 

If p G L{{X}}[y] is any polynomial in Y with coefficients in L{{X}}, the valuation of p at 
7 G L{{X}} is defined to be v.y{p) := v{p{x,'y)). 

4.2. Conjugate Puiseux Series. Two Puiseux series in L{{X}} are called conjugate if they 
are conjugate as field elements over K{{X)). 

4.3. Rational Part. Let 7 = + 02 ^*^ + • • • + akX^>= + ak+iX^’^+^ + • • • G Vx, with 

0 < < ^2 < • • • • Let A: > 0 be such that OjXL g K[X] ioi 1 < i < k and 

Then we call aqX^^ + • • • + a^X^'^ the rational part of 7 , and ak+iX^'=+^ its first non-rational 
term. 
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4.4. Characteristic Exponents. For 7 G Vx, let m G N be minimal with 7 G and 

write 7 = X]j>Q with coefficients bi ^ L. If m = 1, there are no characteristic exponents. 

If m > 2, the characteristic exponents of 7 are defined inductively by 

ei := min{i | 7 ^ 0 and m f i}, 

Cy := min{i | 6 j 7 ^ 0 , gcd(ei,..., e^-i) f z} for > 1 . 

Then ei < 62 <■■■ ■ In fact, there are only finitely many e,^, and these are coprime. 

Example 4.1. If 7 = the common denominator is m = 8 . 

Writing 7 = 2X^1^^X^!^-hX^'^ 1^ , we see that the characteristic exponents are e\ = 4, 
62 = 6 , and 63 = 17. 

4.5. Puiseux Expansions. In what follows, we consider a monic polynomial g G iF[[X]][y] of 
degree m in Y. By the Newton-Puiseux theorem, g has m roots 71 ,... , 7 ^ G L{{X}}: 

5 = (y_7i)...(y_^^)eK[[X]][y]. 

The monic assumption guarantees that each root 7 * is, in particular, integral over T[[^]] and, 
thus, contained in some C Vx- That is, the terms of 7 ^ have non-negative exponents 

only. 

The roots 7 * are called the Puiseux expansions of g {at X = 0). 

4.6. The Newton-Puiseux Algorithm. The Puiseux expansions of g can be computed re¬ 

cursively up to any given order using the Newton-Puiseux algorithm (see, for example, [9]). 
Essentially, to get a solution aiX^^ -|- 02^*2 -|- ... of g{X,^{X)) = 0, with ti < t 2 < ■ ■ ■, the 
algorithm proceeds as follows: Starting from g^^^ = g and = K{{X)), we commence the ith 
step of the algorithm by looking at a polynomial G We then choose one face 

A of the Newton polygon of such that all the other points of the polygon lie on or above 

the line containing the face. Let < 7 ^ be the sum of terms of involving the monomials of 

gr(*-i) on A. That is, if — ^ is the slope of A, then is the sum of terms of of lowest 

( 1 , ^)-weighted degree. We write di for this degree. Choose an irreducible factor of 77 ^ over 

Ai* and a root qi of that factor. Note that Qi is of type qt = CiX^i , where Ci is a root of 
the polynomial g'^^^\l,Y). Now, let = K^'^~^\qi) and set g^^^ = ■^g^^~^\X,qi ■ (1 -|-T)). 
Then the Tth term of the expansion to be constructed is a*A** = Qi''' Qi- It is clear from this 
construction that different conjugacy classes of expansions arise from different choices for the 
faces and irreducible factors of g^^ over respectively. 

Example 4.2. The eight Puiseux expansions of the polynomial 

g = Y^ + {-AX^ + 4A^)y^ + {AX^ - AX^ - 10X^)Y^ + (AX^ - 6 X^)Y^ 

+ ( 6 A® - 8 A®)y^ ( 8 A® - AX^)Y^ + (AX^ -h AX^°)Y'^ + AX^^Y + G Q[A, A] 

are conjugate over Q((A)); their singular parts are of type 

Qi + qiq 2 + qmqs, 

where the qi satisfy 

ql+X^ = 0 , ql + 2 x ^1 = 0 , and qj + j^qi = 0 . 

To see this, note that the Newton polygon of g^^^ = g has only one face Aq, leading to < 7 ® = 
(A^ -|- and the extension 

Ao = Q((A)) c Ai = Ao[iAi]. 

In the next step, g^^'^ has only one face Ai, yielding 

92=4 ( 25 '^ + !)' 
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and 

Ki C K2 = Ko[ixl,{l-i)X^. 
Finally, also has only one face A 2 , which corresponds to 


and the extension 

K 2 CK 3 





Ko[ixl,{l-i)XX{l + i)X^ = 


Ko[i,Xi]. 


4.7. Regularity Index and Singular Part. If 7 = + 02 X ^2 -|-... is a Puiseux expansion 

of g, with 0 < < t 2 < • ■ • and no ai zero, we define the regularity index of 7 (with respect to 

g) to be the least exponent such that no other Puiseux expansion of g has the same initial 
part aiX^^ + • • • + akX^^. This initial part is, then, called the singular part of 7 (with respect 
to g). 


4.8. Maximal Integrality Exponents. Let T = { 71 ,..., 7 ^} be the set of Puiseux expansions 
of g. The valuation of p G L{{X}}[y] at g is defined to be Vg{p) = v-y^ip) ■ Note that 

if p is monic of degree d in Y, where 1 < d < m — 1, and 


p={Y-gy{X))---{Y-gd{X)) 
is the factorization of p in L{{X}}[y], then 

d 

Vgip) = min ^^( 7 * - gj). 

l<t<m 


i=i 


Lemma 4.3. Let g G iL[[X]] [Y] be monic of degree m in Y, with Puiseux expansions 71 ,..., 7 ^- 
Fix an integer d with 1 < d < m — 1. If A C {I,, m} is a subset of eardinality d, set 


lnt(>4) = min V v{-fi - 7 ^) 

\ —■ 

Choose a subset .4, C {1,... ,m} of cardinality d such that Int(.4) is maximal among all Int(^) 
as above, and set p = Wj^^{Y — 7 j) G 'Px\y]- Then Vgip) = Int(^), and this number is the 
maximal valuation Vg{q), for q G L{{X}}[y] monic of degree d in Y. 



Proof. That Vg{p) = Int(.4.) is clear from the definitions. That this number is the maximum 
valuation Vg{q) as claimed follows as in the proof of [21, Theorem 5.1], where the case d = m — 1 
is treated. □ 


In the situation of the lemma, with T = { 71 ,..., 7 m}, we write 

o{r,d) = Vg{p). 

Then 

o(T, 1) < • • • < o(T, m — 1) 

by construction. 

In case d = m — 1, we abbreviate 

Intj = Int({l,..., z - 1 , z + 1 ,..., n}) = ^ z;( 7 i - 7 ^). 

Example 4.4. Let g = (y^ + 2X^) + y^ G Q[y, y]. The Puiseux expansions of g are 

71 = aiy3/2 + + ... , 

72 = a2X^I^ +X^ + ... , 

73 = — 1 — 2 X^ + ... , 
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where ai,a 2 are the roots of + 2. Then Inti = 3/2 + 0 = 3/2, Int 2 = 3/2 + 0 = 3/2, and 
Inta = 0 + 0 = 0, so that both i = 1 and i = 2 maximize the valuation. Taking i = 1, we get 
P= (T- 72 )( 1 ^- 73 ) and o(r, 2) = 3/2. 

Example 4.5. Let g = {Y^ + X‘^){Y‘^ — X^) + T® G Q[X, Y], The Puiseux expansions of g are 

71 = aiX2/3 + ..., ^4 = x3/2 + ..., 

72 = 02^2/3 + . . . , ^5 = -^3/2 + . . . , 

73 = 03^2/3 + . . . ^ 76 = 1 + . . . , 

where the a* are the roots of X3+l = 0. Then Inti = Int 2 = Inta = 2/3+2/3+2/3+2/3+0 = 8/3, 
Int 4 = Ints = 3/2 + 2/3 + 2/3 + 2/3 + 0 = 7/2, and Inte = 0. We conclude that o{r, 5) = 7/2. 

Note that for R any one of the rings K[X], iL[[X]], K{{X)) or Vx, we have R[Y] C 
L{{X}}[y], hence the definition of Vg{p) also applies for p G R[Y]. 

Lemma 4.6. Let g G iL[[X]][y] C L{{X}}[y] be monic of degree m in Y, let 1 < d < m — 1, 
and let R be one of the rings K[X], K[[X]], K{{X)), Vx, or L{{X}}. The maximal valuation 
’OgiQ), Q £ monic of degree d in Y, is independent of the choice of R from among this list. 

Proof. By Lemma 4.3, there is a polynomial p = W^^^iY — 'jj) with |.4| = d which maximizes 
the valuation in case R = L{{X}}, but is contained in Vx[Y]. Choose an integer m such that 
p G L((X^/™'))[y]. By truncating the coefficients of jj to degree Vg{p), we get a polynomial 
p G L(X^/™’)[y]. In fact, for the field K C L obtained by adjoining the coefficients of p to 
K, we have p G iL[X^/™'][y]. Since p is monic in Y, by applying the trace map of 
over K{X) to p and dividing by the integer lead coefficient of the resulting polynomial, we get 
a monic polynomial 0 / p G iL[X][y] of degree d. By construction, p satisfies Vg{p) > Vg{p). 
Note that X[X] is included in all the rings R in the above list, hence the reverse inequalities are 
trivial, and the result follows. 

□ 

Remark 4.7. With notation as above, note that for any representative in K{{X))\Y] of an 
element of iL[X]^x)[^]) the valuation at g is the same. Hence we can also define Vg{p) for 
p G [y] and the maximal valuation Vg{q) with q G K[X](^x)\Y], is the same as the 

maximal valuation with q G iL[X][y], for monic polynomials of the same degree. 

The reason for considering the valuations Vg{p) is that they are directly related to integrality. 
Suppose 

A = K[x,y]=K[X,Y]/{f) 

is the coordinate ring of an irreducible plane curve C of degree n with assumptions and notation 
as in the introduction. 

Definition 4.8. Let R = X[X], K[X\i^x) oi' -^[[^]] B = R[Y]/{f). If G R[Y] is monic in 

y of degree 1 < i < n — 1, and e is the maximal integer such that is integral over B, we 

call eij(g) := e the integrality exponent of q with respect to f and R. 

Lemma 4.9. Let p G K[X](^x) [h"] be monic in Y of degree 1 < i < n — 1. Then 

eK[x]^^^{p) = [vf{p)\- 

Proof. By [19, Theorem 3.2.6], p(x,i/)/x® is integral over K[X]i^x)[^\l if) iff v^iplX^) > 0 for 
every Puiseux expansion 7 of /. Since Vf{p) is defined to be the minimum of the respective 
Vj{p), the result follows. □ 

Lemma 4.10. Let p G iL[[X]][y] be monic in Y of degree 1 < i < n — Then 

ok[[xw{p) = [vf{p)\. 
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Proof. Again, by [19, Theorem 3.2.6], p{x,y)fx^ is integral over A[[A]][y]/ (/) iff v^^p/X^^) > 0 
for every Puiseux expansion 7 of /. □ 

Lemma 4.11. Let p G K[X, Y] he monic in Y of degree 1 < i < n — 1. Then 

eK[x]ip) = [vfip)\. 

Proof. We have to show that p{x,y)/x^ is integral over A iff v.y{pfX^) > 0 for every Puiseux 
expansion 7 of / (see also [21, Section 2.4]). If p{x,y)/x^ is integral over A, then also over 
K[X](^x)\X]/ {f)i hence v^{plX^) > 0 for every Puiseux expansion 7 of /. For the converse, 
note that i;.y(p/A®) > 0 for all Puiseux expansions 7 of / at x / 0. Hence, p{x,y)/x^ G A by 
[19, Theorem 3.2.6]. □ 

Definition 4.12. For q G K[X](^x)\X\ or g G A[[X]][T], we define 

e{q) = [vfiq)\. 

Note that, by Lemmata 4.9, 4.10 and 4.11, we have e{q) = eji{q) as long as g or a representative 
of q is in R[Y]. 

Definition 4.13. Let 0 < i < n be an integer. Taking Lemma 4.6 into account, the number 

e* := max{eR(g) | q G R[Y] monic in Y, degq = i} 

is independent of the choice of R among K[X], K[X](^x) and A[[A]].We call the maximal 
integrality exponent with respect to f in degree i. 

Definition 4.14. We call 

E{f) = en-i = Lo(r,n-l)J 
the maximal integrality exponent of f. 

5. Normalization of Plane Curves via Localization and Completion: 

Decomposing into Branches 

From now on, 

71 = K[C] = K[x, y] = K[X, Y]/{f{X, Y)) 

will be the coordinate ring of an irreducible plane curve C with assumptions as in the in¬ 
troduction. In particular, / is assumed to be monic in Y. We focus on the case where 
P = {X,Y) G Sing(A). Applying the Weierstrass preparation theorem, we get a unique factor¬ 
ization 

( 1 ) f = fofl---fr, 

where /o G A'[[A]][y] is a unit in A"[[y, T]] and fi,..., fr are irreducible Weierstrass polynomials 
in A[[A]][y] (see, for example, [9]). We write 

mi = degy(/i), i = 0 ,...,r 

and refer to /i,as the branches of / at P. 

In this section, we will study first integral bases for the branches of the singularity P from 
a theoretical point of view and explain how these can be combined to give an integral basis for 
K[[X]][Y]/{h---fr). 

In what follows, we consider a monic polynomial g G A'[[A]][y] of degree m in y and write 

B = K[[X]][Y]/{g). 

(We do not assume to be a Weierstrass polynomial, although in this section we apply the 
results to that case.) By abuse of notation, x, y will also denote the residue classes of X, Y 
modulo g. Applying [19, Theorem 3.3.4] to the PID 

K[[x]] C K{{x)) C K{{x))[y] 

we see that H is a free A'[[x]]-module of rank m. 
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Remark 5.1. In particular, i^[[X]][y]/(/i • • • fr) and its normalization are a free -modules 
of rank mi m,.. 

Definition 5.2. With notation as above, we refer to any integral basis of B over iir[[x]] as an 
integral basis for g. 

Lemma 5.3. With notation as above, there exist polynomials pi £ iir[X][y], 0<i<m — 1 of 
degree i in Y with leading coefficients x**, U G Z>o, and Ci £ Z>o, 1 < i < m — 1, such that 

^ Pi{x,y) pm-i{x,y) ] 

is an integral basis for g. Furthermore, ifq £ li'[[X]][y] is any polynomial of degree 1 < k < m—1 
in Y with leading coefficient x^, and e is an integer such that is integral over -fC[[x]], then 

e < Ck- In particular, the e* depend only on g and satisfy 0 = eo < ei < • • • < Cm-i- 

Proof. Each square matrix with entries in ii"[[x]] of maximal rank has a uniquely determined 
upper triangular Hermite normal form {pij), where the diagonal elements are of type pu = 
and where the pij, j > i, are polynomials in K[x] of degree < r'j (see [12]). Hence, given any 
integral basis for g where the denominators are powers of x, we can first reduce the numerators 
modulo the monic polynomial g to get elements of y-degree at most m — 1. Then taking the 
largest power of x in the denominators as common denominator, we construct the matrix of 
coefficients of the numerators and apply unimodular row operations as in Remark 1.2 to get 
after cancellation an integral basis 1 = 5o, 6i,... 6m-ij where each bi is of the form ; with 

ei,ti £ Z>o {ti < Vi), Pi £ ib[y][y] polynomial of degree i m Y with leading coefficient x*k 
This shows the first statement of the lemma. The second statement and, thus, the uniqueness 
result follows by expressing as a ib[[x]]-linear combination of the bi, 0 < i < k. To see 

that ej-i < ei, for each i, consider q = Y ■ pi-i. □ 

Remark 5.4. We refer to [12] for the computational aspects of the lemma up to any desired 
precision (that is, up to which power in X the coefficients are developed). We note that in our 
case, since the starting point is an integral basis, we know that l,y,..., can be expressed 
as ib[[y]]-linear combinations of the elements in the basis, hence the exponent of the common 
denominator of the input basis gives an a-priori bound for the maximum precision needed. 

Remark 5.5. We will show after the next proposition that, in fact, any integral basis in the 
shape of the lemma is guaranteed to have monic numerators, that is ti = 0 for 0 < z < m — 1. 

In Section 7.5, we will present a practical method for finding integral bases as in Lemma 5.3. 
The starting point for this is the following proposition. 



Proposition 5.6. With notation as above, for 1 < i < m — I, let monic polynomials pi £ 
ib[[X]][y] of degree i in Y be given, and let Ci be the maximal integrality exponent with respect 
to g in degree i. Then 

^ Pi{x,y) pm-i{x,y) 

xe(pi) ’ ■ ■ ■ ’ 

is an integral basis for g iff e{pi) = e* for 1 < i < m — 1 . 


Proof. If the given elements form an integral basis for g, then necessarily e{pi) < Ci for each i. 
Suppose that e{pi) < Ci for some i, and choose an element q £ ib[[X]][y] which is monic in Y 

which is impossible 


of degree i and satisfies e{q) = Cj. Then -^ £ (l 


pi{^,y) 


Pi{^,y) 


K[[x]] 


since the exponent of each denominator on the right hand side is smaller than e, by the third 
part of Lemma 5.3. 

For the converse, suppose that e{pi) = e* for each i. Set B' = /l, • • •, 

Then B <Z B' <Z B, and we have to show that B' = B. That is, given a polynomial q £ Rr[[y]][y] 
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of degree 0 < i < m - 
G B[, where B'- 


1 in y such that SfeMl 
_ /i pi{^,y) Pt{x,y ) ' 


G B for some integer e > 0, we have to show that 
'k[[x\\ 


We do induction on i. There is nothing to show in case i = 0. If i > 1, let c be the leading 
coefficient of g G iir[[X]][y]. We can assume that c = x*, t G Z>o, since any other factor is 
invertible in itr[[W]]. Write g as a product q = x^q, with q G K{{X))[Y] monic in Y. By Lemma 
4.6 and the definition of pi we have Vg{^ < e{pi), hence 


e < e{q) < t + Vg{^ < t + e{pi) = e{x^pi). 


This implies that ^ P^(^^y) g hence it is in B'-. Since degy(g — x^pi) < 

X Pii^^y) ^ ]^ \jy induction hypothesis, we get ^ G B'-_^ C 

G B' as claimed. 


i and _ 

B[. Therefore 

□ 


Remark 5.7. Together with Lemma 4.6, the last proposition proves the existence of an integral 
basis |l = poj , • • •, where the pi G iL[X][y] are monic of degree i in T, 0 < i < 

m — 1. 

Proposition 5.8. With notation as above, i/ |l = po, , • • ■, } is any integral 

basis of g in the shape of Lemma 5.3, the polynomials pi in the integral basis are monic in Y. 


Proof. Suppose for pk we have tk > 0. Since no cancellation is possible, there must be some 
coefficient of pk that is not multiple of x. We now take an integral basis B for g as in the last 
remark, and call Bk the elements in B of degree at most k. We can express as a 

iL[[W]]-linear combination of the elements in Bk (no elements in B of larger degree can be used). 
But the largest power of x in the denominators of Bk is e^, hence we get a contradiction. □ 


We now return to the branches /i,..., /^ at P of our given polynomial / and apply the above 
to the product g = fi ■ ■ ■ fr- 

Proposition 5.9 (Splitting of Normalization). Let fi,..., fr be the branches of f at P as in 
Equation (1). For each i, set hi = 0^=1 Then the fi and hi are eoprime in K{{X))[Y], 

so that there are elements ai,bi G iL[[X]][y] and integers Cj G N such that 

a-ifi + bihi = X^\ for i = 1... ,r. 

Furthermore, the normalization of K\fX\\\Y]/{fi ■ ■ ■ fr) splits as 

r 

K[[X]][Y]/{h •••/.) = 0 K[[X]][Y]/{f,), 

i=l 

and the splitting is given by 

{ti mod fi,...,tr mod fr) ^ mod /i • • • /r- 

i=l 

Proof. Clear by the Chinese remainder theorem and its proof. See [9, Theorem 1.5.20]. □ 

Given an integral basis for each branch, we can make the splitting of normalization explicit: 


Corollary 5.10. With notation as in Proposition 5.9, for i = 1,..., r, let 


l=Po 


(0 

» Pi 




(0 

Prn,-1 


yep 

represent an integral basis as in Lemma 5.3 for fi, and set 




bihi bihipl 




bihip^^. 




(0 


X 


I „(i) 
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or each i. Then U ... U is an integral basis for /i • • • /r • 

Proof. Immediate from Proposition 5.9. □ 

Finally, we note that we can apply the construction in Lemma 5.3 to the elements in B^^'^ U 
... U to get an integral basis in the shape 

_ Pi{x,y) pm-i{x,y) 

with polynomials pi G iF[X][y] of degree i in y and non-negative integers ej,l<z<m — 1, 
where the polynomials pi are monic by Proposition 5.8. 

Using some of the tools developed in the subsequent sections, we illustrate the corollary by 
an example: 

Example 5.11. Let / = (y^ -|- X^)(y^ — X^) + y® be as in Example 4.5 and let A = 
K[X,Y ]/(/) = K[x,y]. In iL[[X]][y], / can be factorized as / = / 0 / 1 / 2 , whose developments 
up to degree 3 in X are /o = E + {-X^ - X^ + 1), fi = Y^ + {X^ + X‘^)Y^ + {-X^)Y + 
and /2 = y^ — X^. 

Following Proposition 5.9, applying the extended GCD algorithm to fi and hi = f^ we get 
the coefficients ai and bi whose developments up to degree 3 in X are ai = —4X^y — 2X^ — 
2X^Y - X'^ + XY -Y - 1 and bi = -dX^Y'^ - 2X^Y - 2X‘^Y'^ - 3X^ - 2X'^Y + - E^ - y 

satisfying ai/i + bihi = X^, with hi = f 2 - 

The rings X[[X]][y]/(/j), i = 1,2, have integral bases |l,y, and {l,|}- By Corollary 

5.10, an integral basis basis for the normalization of X[[X]][y]/(/i/ 2 ) is given by B^^^ U B^'^f 
where 

bif2 bif2y bif2y‘^\ ^( 2 ) ^ «i/ii/ ] 

\ ’ x^ J 


BA) = 


x^ 


We can now apply the construction from Lemma 5.3. Since the maximum power of x appearing 
in the denominators of B^^^ U B^'^^ is X3, we can truncate all the coefficients appearing in the 
computation to degree 3 in x. We obtain the integral basis 

y2 y3 y^ + x'^y 

; I/; ; o ? q 

X X^ X'^ 

for the normalization of X[[X]][y]/(/i/2). 

Note that in this example, the maximum power of X appearing in the denominators was not 
known a priori, but in a practical algorithm it is required for computing the developments of 
the factors /i, /2 and the coefficients ai, bi. We address this problem in Proposition 7.10. 


6. Normalization of Plane Curves via Localization and Completion: Local 

Contributions 

In this section, we will keep the notation and assumptions of the previous section. For 
simplicity, we will assume that the origin is the only singularity at X = 0, which can always be 
achieved by a linear coordinate change. Alternatively, it is easy to extend our algorithms for 
the case of the presence of more than one singularity at X = 0. 

Using Puiseux series, we will show how to pass from an integral basis for the normaliza¬ 
tion of X[[X]][y]/(/i • • • fm) as in Corollary 5.10 to an integral basis for the normalization of 
X[[X]][y]/(/). If the elements of that basis are polynomials, we will show that this is already 
an integral basis for the normalization of K[X](^x)[^]/if)■ Moreover, since we are assuming 
that the origin is the only singularity at X = 0, it will also be an integral basis for the local 
contribution to the normalization at P = (X, E). 

Proposition 6.1. Let f = fog be a factorization of f with fQ,g G X[[X]][y], /o a unit in 
X[[X]][y] and g a Weierstrass polynomial, let A = X[[X]][y]/ (/) and P = (X, E). Write 
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n = degy f, do = degy /o and m = degy g. Suppose that P is the only singularity of f at 
X = 0. Let 


_ ^ Pm-l 

be an integral basis for iC[[X]][y]/( 5 r) with polynomials pi G i^[X][y] monic in Y of degree i. 
Let fo £ K[X][Y] be a (monic) polynomial with 


fo = /omodX®’"-! 


Then 


\/oPo, 


foPi 


is an integral basis for the normalization of K^X\\\Y]/{f). 


foPm—l 


Proof. We first show that the elements foPijx^L 1 < i < m — 1, are integral over A. Suppose 
that 7 is a Puiseux expansion of /. If 7 ( 0 ) = 0 then v.y{foPi) > v^{pi) > Vg{pi) > Ci. If 7 ( 0 ) 7 ^ 0 
then v^ifoPi) > v^ifo) > e^-i > e* by definition of /q. 

To apply Proposition 5.6, we have to show for i = 0,..., m — 1 that Cj is maximal among all 
e{q) where q G iir[[X]][y] monic in Y of degy(g) = i + do, and that e{q) = 0 for all q G i^[[y]][y] 
with degy(g) < do- Suppose q{x,y)/x^, e > 0 is integral over A where q G iir[[y]][y] is monic 
in y of degree s < degy(/). Let 

q = qoh 

be a factorization of q with qo,h G iL[[X]][y], qo(0,0) not 0 and h a Weierstrass polynomial. 
Write 7 / 1 ,..., 7 s for the Puiseux expansions of q. 

The assumption that the origin is the only singularity at X = 0 implies that all the Puiseux 
expansions of fo have pairwise different non-zero constant term. Suppose that 7 ^)^... ^ 7 (‘^o) 
are the Puiseux expansions of fo and ,..., 7 ("^) are the expansions of g. Then 7 ^*^ = 

Og ^ '+..., 1 < i < do) with Og ^ / 0 for all 1 < i < do and Og ^ for 1 < i < j < do. 

Moreover, -l-... for i > do- Assume that for some 1 < i < do there is no expansion 

(i) 

dji 1 ^ i ^ 'S) with initial term ay. Then 

S 

Vf{q) = mm v ii){q) = mm y'u(7i - pj) = 0, 

l<2<n ' l<i<n 

which contradicts the hypotheses. 

(i) 

Hence, for s < do the maximal integrality exponent is 0. Moreover, if s > do then all Og' 
for 1 < i < do have to appear as initial terms of some gj, 1 < j < s, and they are pairwise 
different. So we can assume that for each 1 < i < do the initial term of 7 * is Og^ Then 
Qo = {y — Pi) ' ■ ■ ■ ' {y — Pdo) ■ where u G Vx\Y], in particular 


do < degy (go). 

For any Puiseux expansion 7 ^*^ of / we have 

( \ I (hi I if7^*Ho) = 0 

( 9 ) = (®) + (A) = I . J ^ 

Vf{q)= min u (i) (g) = min min V (i){qo), min V (i){h)} 
l<i<n ' ' do<i<n ' J 

= mm{vfg(qo),Vg(h)} < Vg(h). 


hence 


This implies that 

e < [vf{q)\ < [vgih)\ < es-deg^igo)- 

For the last inequality, we use that es_degy(qo) maximal among all [vg{p)\ with p G iF[[X]][y]. 
Since do < degy ( 70 ), we have 

6 y C5_degy(gQ) y 
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which proves our claim. □ 

Remark 6.2. For the last proposition, we do not use the assumption pi G 0 < i < 

m — 1, and we do not need to truncate /q. However, by doing this we obtain an integral basis 
with numerators in iF[X][y], which is then, by faithfully flatness, also an integral basis for the 
normalization of K[X](^x)[^]/if)■ Moreover, if P = {X,Y) is the only singularity at X = 0, it 
is also a X[X]-basis for the local contribution to the normalization of K[X, Y]/{f) at the origin, 
as we show next. 


We need first the following version of Lemma 5.6. 

Lemma 6.3. Let f ^ K[X,Y] monic of Y-degree n Let pi,... ,pn-i € K[X,Y] be polynomials 
which are monic in Y of degree i such that for all 1 < i < n — 1, e{pi) is the maximal integrality 
exponent with respect to f in degree i. Then 

^ Pi{x,y) pn-i{x,y) 

xcipi) ’ ■ ■ ■ ’ j,e{pn-l) 

form a K[x]-module basis of a local contribution to the normalization at any prime Q with 
{x) C Q. If P = {X,Y) is the only singularity at X = 0, then it is a K[X]-module basis for the 
minimal local contribution at P. 


K[x] 


as we 


Proof. Suppose that e{pi) is the maximal integrality exponent with respect to / in degree i. By 
Corollary 4.6 

e{pi)> Vf{p) 

for all p G K{{X))[Y] of degree i. 

Denote by A' the iF[x]-module generated by 1, • ■ •, • By assumption, W C A. 

We prove that A' the minimal local contribution at P. 

First, given a polynomial q{X, Y) G K[X, y] of degree Q < i < n — 1 vnY such that G A 

for some integer e > 0, we show that G A[, where A[ = ^1, ■ • •, 

did in the proof of Lemma 5.6. 

We do induction on i, the claim being trivial for i = 0. For i > 1, dividing q by the lead 
coefficient with respect to Y write q{X,Y) = c{X)q'{X,Y) with c{X) G K[X] and q'{X,Y) G 
K{{X))[Y], monic in Y. By the above remark and the definition of pt we have Vf{q') < e{pi), 
hence 

e < e{q) < e(c) + Vf{q') < e(c) + e{pi) = e{cpi). 

(Here the last equality holds, since c{X) has integer valuation and only depends on X). This 
implies that Moreover, writing c'(x) = we obtain that 

c{x)pi{x,y) , sPi{x,y) , 

-:- = c [x) —^— 77 -^ G A'. 


Since degy(g — cpi) < i and — 

q{x,y) 


x^ ' ' 

c{x)pi{x,y) 


G A, by the induction hypothesis, we get 
c{x)pi{x,y) 


G C A-. 


x^ 


Hence, G A'^ as claimed. 


We now pass to the localization. For 

D = K[X]^x)[Y\/{f) and D'= 


= 1 


Pi{x,y) 


x' 


e(Pi) 


Pn-i{x,y) \ 

XpPn — l) j 


K[x\^^) 


we have D C D' C D. Let G D be an arbitrary element of D with polynomials q G K[X, y] 


d{x) 


and d G K[X]. Write d{x) = x® • d'{x) with a unit d'{x) G K[x](^j.y Then also 
there exists an h{X) G K[X] such that h{x) G K[x](^^'^ is a unit and h{x 
above argument, /i(x)^^||^ G A C D', so G D', hence, D' = D. 


g{^,y) 


G D, hence 




G A. By the 
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To localize N at primes of A, we first prove that Pl is a commutative ring with 1. This 
amounts to show that N is closed under multiplication. Any product of elements can be written 
as 

g(a:,?/) _ 4{x,v) ^ q''{x,y) 

with q" G K[X,Y] and degy(g") < n. Since the product is in A, by the above argument, the 
product is in A! . Since A' is a ring it is also an A-module. 

We have 


^'q = Dq = Dq = Aq 

for all Q G Spec A with (x) C Q. Moreover, A'q = Aq for all Q with (x) ^ Q, since the 
denominators of the generators of A' are in (x). If P = {X, Y) is the only singularity at A = 0, 
then Aq = Aq for all Q G Spec A with (x) C Q, Q P, hence it is a minimal contribution. □ 


Corollary 6.4. Let 


B = <l=po, 


Plix,y) Pn-lix,y) 


be an integral basis for K[[X]][Y]/(f), where pi are polynomials. Then it is also an integral basis 
for the normalization of K[X](^x)\^]/{f) a K[X]-module basis for the local contribution at 
any prime Q with (x) C Q. If P = {X,Y) is the only singularity at X = 0, then it is a 
K[X]-module basis for the minimal local contribution at P. 


Proof. Clear by the previous lemma and Lemma 4.6. 


□ 


7. Normalization of Plane Curves via Localization and Completion: The 

Algorithmic Point of View 

Let A = K[x, y] = K[X, Y]/{f{X, Y)) be as before. In this section, we show how to compute 
a local contribution to A at each prime ideal P G Sing(A) via Puiseux expansions, Hensel’s 
lemma, and Proposition 5.9. The normalization A itself and an integral basis for A over K[x], 
respectively, are then obtained along the lines of Proposition 3.1 and Remark 1.3. 

We start with a sketch of the algorithm. 

7.1. Summary of the Algorithm. From a theoretical point of view, the algorithm involves 
the following steps: 

(1) If the prime ideal P G Sing(A) corresponds to a (single) AT-rational singularity, translate 
the singularity to the origin. If P corresponds to a set of conjugate singularities over 
K, extend the base field K as needed, and translate one of the singularities to the 
origin. In any case, apply a linear transformation so that the translated singularity is the 
only singularity at A = 0. 

For the singularity at the origin, do (to simplify the presentation, we will still write / for the 
transformed equation of our curve): 

(2) Determine the maximum integrality exponent E{f) as described in Section 4.8. 

(3) Determine integers Cj, 1 < z < r, as in Proposition 5.9 and the factorization / = ni=o 
of / into branches, developing each /j up to degree E{f) + c* in A. Here, make use of 
Hensel’s lemma as described in Sections 7.3 and 7.4. 

(4) Compute the Bezout coefficients bi, 1 < i < r, from Proposition 5.9 up to order E{f) + c*. 

(5) For each branch /i, 1 < ^ use the algorithm from Section 7.5 to compute integral 
bases for the branches A[[A]][y]/(/i) as in Lemma 5.3 up to order E{f) + a in A. 

( 6 ) Construct the generating sets 1 < z < r, as in Corollary 5.10, and apply the construc¬ 
tion in Lemma 5.3 to compute an integral basis for A[[A]][y]/(/i • • • fr) with numerators 
in A[A,y]. 

(7) Compute an integral basis for A[[A]] [Y]/{f) using Algorithm 10, which is based on Propo¬ 
sition 6.1. Since the numerators in the output are elements in A[A][y], by Corollary 6.4, 
this is the local contribution to the normalization of A at the origin. 
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(8) Apply the inverse translation to the elements of the local contribution to restore the 
singularity to the original position. 

( 9 ) If P corresponds to a set of conjugate singularities, then use Remark 7.17 to modify 
the numerators and denominators of the local contribution in order to obtain the local 
contribution to A at P over the original field. 

From a practical point of view, we face the problem that, in the approach outlined above, 
we need to determine the Cj a priori. Moreover, the computation of the Bezout coefficients bi 
via the extended Euclidean algorithm is very time consuming. To remedy these issues, relying 
on Proposition 7.10 below, we will replace the bi and Cj in Steps 4 and 3 by easier to construct 
polynomials Pi £ K[X,Y] and appropriate vanishing orders, respectively. 

We refer to the following sections for more details. 

7 . 2 . Puiseux Expansions. The factors fi appearing in the decomposition 

f = fog = fofi ■■■ fr 

of / into its branches and a unit /o as in Equation ( 1 ) in Section 5 correspond to complete sets 
of conjugate Puiseux expansions. Developed up to a given degree, the fi may hence be found by 
computing the expansions and their respective products. Since this is computationally involved, 
we propose a different approach which, via Hensel’s lemma, makes considerably less use of the 
Newton-Puiseux algorithm. In describing the new approach, we use the notation below. 

We partition the set of all Puiseux expansions of / into Puiseux blocks. A Puiseux block 
represented by an expansion 7 with 7 (0) = 0 is obtained by collecting all expansions whose 
rational part agrees with that of 7 and whose first non-rational term is conjugate to that of 7 
over K{{X)). A Puiseux segment is defined as the union of all blocks having the same initial 
exponent. That is, we have one Puiseux segment for each face of the Newton polygon of /. In 
addition, all Puiseux expansions 70!/ with 7 ( 0 ) / 0 are grouped together to a single Puiseux 
block of an extra Puiseux segment. In this way, the Puiseux expansions of / are divided into 
Puiseux segments, each segment consists of Puiseux blocks, and each block is the union of classes 
of conjugate expansions. 

Example 7.1. Suppose that the Puiseux expansions of the given polynomial / are 


71 = 1 + ^2 + . . . , 

76 = X + 6 iX 3/2 + x 3 + 

72 = -1 + 3 A + ..., 

77 = X + 62^3/2 + X 3 + 

73 = aiA 3/2 + 2X2 + ..., 

78 = X + 5 iX 3/2 + X^ + 

74 = 02X3/2 + 2X2 + ..., 

79 = X + 62x3/2 + X^ + 

75 = X + 3 X 2 + ..., 



where {73,74}, {76,77} and {73,79} are pairs of conjugate Puiseux series. Then {71,72} is the 
segment of expansions 7 with 7 ( 0 ) 7^ 0. Another segment is {73,74} (which consists of one 
block containing a single class of conjugate expansions). All the other expansions form a single 
segment, consisting of the blocks { 75 } and {76,77,78,79}- The last block contains two classes 
of conjugate expansions, namely {76,77} and {73,79}. 

7 . 3 . HensePs Lemma. We begin by recalling the statement of the lemma: 

Lemma 7.2. Let F £ Ar[[A]][y] be a monic polynomial in Y, and assume that F{0,Y) = goho, 
with monic polynomials go, ho £ K[Y] such that {go, ho) = A"[y]- Then there exist unique monic 
polynomials G,H £ Ar[[A]][y] such that 

(1) F = GH, 

( 2 ) G{ 0 ,Y)=go, HiO,Y) = ho. 

In fact, for each m G N, there exist unique gm,hm £ K[X,Y] of X -degree < m such that 

(3) F = gmhm in {K[[X]y{X^+^))[Y], 
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(4) gm^9i, hm^hi in {K[[X]]/{X^+^))[Y], i = 0,... ,m - 1. 

Proof. See, for example, [1]. □ 

Conditions (3) and (4) imply that the polynomials gm and can be computed inductively 
along the X-degree, solving for each m a system of n linear equations in n variables, where n is 
the y-degree of F: For each 0<z<n — 1, we get an equation by comparing the coefficients of 
Xmy^ in F and in gmhm- For further reference in this paper, we state the resulting procedure 
as Algorithm 1, omitting the actual computation steps. 


Algorithm 1 HenselLift 

Input: F £ K[X,Y] monic in Y] go, ho G K[Y] monic with F{0,Y) = goho, {go, ho) = A[y]; 

den. 

Output: g,h £ K[X,Y] of A-degree < d, with g{0,Y) = go, h{0,Y) = ho, and F = gh 
mod X^^^. 


When applying HenselLift as indicated in Section 7.2, we hrst address the Puiseux segment 
consisting of all Puiseux expansions 7 of / with 7 (0) / 0. That is, we decompose / as / = fog 
as in ( 1 ), separating the unit fo from the component g vanishing at the origin (we develop fo 
and <7 up to a certain order). This is summarized in Algorithm 2. 


Algorithm 2 SeparateUnit 

Input: / £ K[X,Y] irreducible and monic in Y, with /(0,0) = 0; d G N. 

Output: fo,g £ iF[[A]][y] as in Equation (1), developed up to A-degree d. 

1 : compute monic go, ho £ K[Y] with ho = Y^ for some k £ N>i, Y f go, and /(O, Y) = goho 
2 : return HenselLif t(/, go; ^0; 


Example 7.3. Let / = (y — A)(y+A)(y + 2A)+y^. Then there are three Puiseux expansions 
satisfying 7 ( 0 ) = 0 and four expansions with 7 ( 0 ) 7 ^ 0 (note that f{0,Y) = y^ + y'^ = y^(l + 
y^)). We write 71 ,... ,74 for the latter expansions and 75 = A + ... ,76 = —A + ... ,77 = —2A + 
... for the expansions vanishing at the origin. We apply SeparateUnit to develop the products 
fo = 7i • • • 74 and 9 = 75767 ? up to degree 2. It calls HenselLift( f, on, ho, 2) with go = I + Y^ 
and ho = Y^. The output is 92 = 5X^Y^ - 2XY^ + y^ + 1 and h 2 = Y^ + 2XY^ - 2X^Y. 

Alternatively, we could compute / = fog by means of the Weierstrass Division Theorem. 
However, the use of Hensel’s Lemma allows for more generality since it does not require to have 
one factor vanishing at the origin. This will be useful for our local version of Hensel’s Lemma. 
It is also useful when the singularity has no A-rational coordinates, as we can modify our 
algorithms to avoid moving the singularity to the origin, which requires the use of an algebraic 
extensions. (For keeping the presentation clear, we do not give the details of this strategy.) 

7.4. A Local Version of HeusePs Lemma. Having computed the decomposition / = fog as 
in the previous section, our next goal is to separate the different Puiseux segments corresponding 
to g. Here, we cannot apply Hensel’s lemma directly since all factors of g vanish at the origin, 
so no matter how we choose go, ho, the condition {go, ho) = K[Y] will not be satisfied (consider, 
for example, the product {Y — 7 i)(y — 72 )(y — 73 ) in Example 4.5). 

To overcome this problem, we transform g as explained in what follows. Write 

71 = uJA*! + a^A*^ ^ 

72 = a? A*? + aiA*i + ..., 


7, = alX^‘ + a|A*i + ... 
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for the Puiseux expansions of g, and suppose for simplicity that t := t\ = mini<j<st5_- Naively, 
we are now tempted to substitute X^Y for Y in g = (y — 71) • • • (y — 7^) G K[[X]][y] and cancel 
out X* in all factors in order to separate the Puiseux segment corresponding to t from the rest. 
However, since this would introduce fractional exponents and thus force us to leave it'[[X]][y], 
we proceed in a slightly different way: Write t = u/v, with u,v £ N>i, and substitute X'^ for X 
and X“y for Y. Then set 

F(x,y) = 5(x^x“y)/x"“ 

= (y - (al + + ...))••• (y - (a^X^f + ...))€ X[[X]][y]. 

Now, F has factors not vanishing at the origin, and these correspond to the Puiseux expansions 
of / forming the Puiseux segment with smallest initial exponent t. Applying Hensel’s lemma, 
reversing the transformation, and iterating the process yields Algorithm 3. 


Algorithm 3 SegmentSplitting 

Input: g G X[X, y] irreducible and monic in Y, with <7(0,0) = 0 and no Puiseux expansions 
vanishing at y = 0 outside the origin; d G N. 

Output: Weierstrass polynomials gi, ■ ■ ■, gi £ X[[X]][y], all developed up to X-degree d, with 
g = gi ••• Oi ss in (1), and each gi corresponding to precisely one Puiseux segment of g as 
outlined above. 

1: let be the different initial exponents of the Puiseux expansions of g (which are 

obtained from the Newton polygon of g) 

2: if £ = 1 then 

3: return {g} 

i- t = u/v = min{ti,..., t^}, with u,v £ N 

5 : s = y-degree of g 

6- F = c/(X’^,X“y)/X*“ 

7: compute Gq, Hq £ K[Y] with F[q = Y'^ for some w £ N>i, Y f Gq, and F(0, y) = GqHq 

8: {G, H} = HenselLift(F, Go, 

9: gi = G(X^I'".YjX'^l'^), h = H'(XV»^,y/X“/^) 

10: return {^ri} U SegmentSplitting(/i) 


See [9, Theorem 5.1.17] for an alternative approach extending the Weierstrass Division The¬ 
orem. 

Example 7.4. Let / = (y2 + 2X3)((y + 2 X ^)2 + X^) + y^. Evaluating / at X = 0, we get 
/(O, y) = (y^ -|- l)y^. Applying SeparateUnit(/, 8 ) gives 

fo = -48X®V + 210X® + 8X’'y - 56X^ - 32X‘*y + 4X'’ + - X’^ + 12X^ - 2X^ - 4X^Y + Y^ - 1, 

g = 48X^Y^ - 46X®y^ + 8X’'y^ + 16X®y + 8X’'y^ - 32X®y^ + 2X® + 4X®y^ 

+ 8X®y^ + ax’" + x®y^ + sx^y + 4x"‘y^ + 2 X^y^ + 4x^y^ + y'‘. 

The Puiseux expansions of g are 

71 = aiX3/2 + aiX^/^ _ ^ , 

72 = 02X3/2 + 02X9/2 - 4X3 + ... , 

73 = -2X2 5^x 3/2 165^x33/2 + . . . ^ 

74 = -2X2 ^ ^ 1662X33/2 + . . . ^ 

where oi, 02 are the roots of X 2 -|- 2 and 61 , 62 are the roots of X 2 -|- 1 . 

The smallest initial exponent t of the expansions is t = u/v = 3/2. We compute 

F{x, Y) = g(x^, x^y) = 48X^^y^ - 8x”y^ + 46X^°y^ + 32X®y'^ + 8X®y^ + 8x’'y^ 

- lex^’y - 4x®y^ + x'*y^ + 2X‘‘ + 4x^y^ + 4xy^ + y^ + sx^ + sxy + 2Y^. 
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Now, F(0,y) = (y^ + 2)y^. Applying Hensel’s lemma to the factors + 2 and Y^, and 
reversing the transformation, we obtain 

gi = -4X^ - 8X^Y + 2X^ + Y^ = (Y^ + 2X^) - 8X^Y - 4A®, 
h = X^ + 4A^ + 4X^Y + y2 = (y + 2X^)2 + 

Next, we wish to separate the different blocks in a given Puiseux segment. Consider first 
blocks inside a segment which have the same initial exponents but whose initial terms are not 
conjugate. In this case, after applying the above transformation and substituting 0 for X, 
F(0, Y) will have different factors that do not vanish at the origin, corresponding to each of 
these blocks. Hence we can still separate these blocks using Hensel’s lemma as before. 

To be able to separate all blocks, it remains to consider the separation of blocks that have the 
same initial rational term (and therefore the same initial exponent). Suppose that gi is a factor of 
/ containing some Puiseux blocks of / all having the same initial terms g = + • • • + a^X ^^, 

a, G A, G N, 1 < i < /c. In this case, we first apply the transformation Y = Yi Y g, and 
compute 5 i(X, Yi) = gi{X, Yi+g). Then gi will contain the same expansions as gi but without 
the initial terms g. We can now proceed as before to separate the blocks. After computing the 
factors corresponding to each block, we replace Ti by T — r/ to get the desired factor. 

Algorithm 4 summarizes this strategy. In Line II of the algorithm, the presence of a power of 
a linear factor implies that the expansions share a common rational part, and hence it is possible 
to further split the factor. 


Algorithm 4 BlockSplitting 

Input: g G K[X, Y] monic of T-degree s > 0 such that all the Puiseux expansions of g are in 
the same Puiseux segment; d G Nq. 

Output: gi,... ,gr G K[X, Y] such that the expansions of each gi are the same as the expansions 
of the i-th Puiseux block of g up to order d in X. 

1: L = 0 

2-. g = the common rational part of all the Puiseux expansions of g 
3: if r/ = 0 then 

4 : t = u/v the initial exponent of the Puiseux expansions of g (which is obtained from the 

Newton polygon of g, and by assumption is the same for all expansions) 

5: g{X,Y)=g{X\X^Y) 

6: F = g/X^^ 

7: Compute Gq, Hq G K[Y] with Gq ^ 1 irreducible or a power of an irreducible polynomial 

and Go, Hq coprime such that F{0, Y) = GqHq. 

8: if FIq 7 ^ I then 

9 : (G, FI) = HenselLif t(A, Gq, Hq, vd) 

10: gi = G(X^/’', y/X“/’'), h = H{X^I^, YjX^/^) 

11: if Go is not a power of a linear factor in Y then 

12: return {gi} U BlockSplitting(/i) 

13: else 

14 : return BlockSplitting(g'i) U BlockSplitting(/i) 

15: else 

16 : return {g'\ 

17: else 

18 : g = g{X,Y + g) 

19 : {gi, ...,gr} = SegmentSplitting(5) 

20: for 1 < i < r do 

21: ,..., /19I = BlockSplitting(g,) 

22: L = LY{hi{X,Y -g),..., hs{X, Y-g)} 

23: return L 
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The ideas from [9, Theorem 5.1.20] can in some cases also be used for our purpose. However, 
the cited theorem is not as general as we require. 

Our final goal is to separate all factors corresponding to different conjugacy classes of expan¬ 
sions. For this, all algorithms known to us require that we work in algebraic field extensions. 
We compute the conjugate Puiseux expansions 71 ,...,% up to the required degree and then 
compute the product (T — 71 ) • • • (T — %). (See Algorithm 10.) This last step is only needed 
when a Puiseux block contains more than one conjugacy class of expansions. 

In Algorithm 5, we sum up the discussion above, arriving at a general splitting algorithm. 

Algorithm 5 Splitting 

Input: / G K[X,Y\ irreducible polynomial, monic of T-degree n; d G Nq. 

Output: /o as in (1) and Weierstrass polynomials G Ar[[A]][y], all developed up to 

A-degree d, with f = fofi ■ ■ ■ fr as in ( 1 ), and each fi,l < i < r corresponding to precisely 
one conjugacy class of Puiseux expansions of / at the origin. 

1: {/ojfl'j = SeparateUnit(/,d) 

2 : {gi,..., gs} = SegmentSplitting( 5 f, d), the factors corresponding to the different Puiseux 
segments of g 

3: for all i = 1,..., s do 

4: compute {hi ,..., hg'} = BlockSplitting( 5 i, d) 

5: for j = 1,..., s' do 

6 : Compute T = { 71 ,..., 7 ^}, the singular part of the expansions of hj 

7: m = number of conjugacy classes in F 

8 : if m > 1 then 

9: for k = 1,... ,m do 

10 : Compute Tfc = { 7 ^, 1 ,.. • , 7 A:,Sfe}) the expansions of the k-th. conjugacy class of T, 

up to order d in A 

11: Pk = {y - 1k,i) ■ ■ ■ {y - lk,sk) developed to degree d 

12 : L = LU{pi,... ,pk} 

13: else 

14: L = LU {hj} 

15: return L. 


7.5. Integral bases for the branches. Let g G A[[A]][y] be an irreducible Weierstrass poly¬ 
nomial of degree m. We show how to algorithmically compute an integral basis for A[[A]] [y]/( 5 f) 
over A[[A]]. That is, we compute polynomials pi, ..., Pm-i £ described in Propo¬ 

sition 5.6 and their corresponding integrality exponents. 

For each d, 0<d<m — 1, we look for a polynomial pd G A [A] [A] of degree d with maximal 
valuation at g. 

Let r be the set of Puiseux expansions of g. Since we are assuming g is irreducible, all the 
expansions of g are conjugate. 

For any d G Nq, 0 < d < m, note that 


o(r,d) = max ng Hiy - v) 
#N=d [ \v&N 

Recall that for a given N C T, we have the formula 


^9 ( 

\t]&N j ' {ri^N J 

To compute o(T, d), 1 < d < m, we do not apply the above formulas but we compute a 
polynomial pd G K\X,Y\ of A-degree d such that Vg{pd) = o{r,d), recursively truncating the 
expansions of g. 
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We consider first the simple case when there exists t € Q such that the conjugated expansions 
7 i,..., 7 m of g agree in the terms of degree lower than t and have conjugate terms G Vx, 
that is 

7i = nix'll + aaX'^^ + ... + ^ ^ ^ _ 

where aj G K and dj G No, 1 < j < /c. To compute the numerator pd of the element of degree d 
in the integral basis, we truncate 7 i to 7 * for 1 < i < d to degree dk and we set 

Pd = {Y-ji)---{Y-jd)€K[X,Y]. 

Lemma 7.5. The polynomial pd defined above has maximal integrality exponents among all 
monic polynomials of degree d in Y. 

Proof. Let = (X — 7 *J ... (X — 'jifi G Vx[Y], 1 < ii < • • • < irf < m, be an element of degree 
d in X of largest valuation at g. We know that we can always take pd in this way. Let i' be an 
index not appearing in {ii,..., id}. We have by construction 

d d 

VgiPd) = v^., {pd) = u( 7 i/ - 7 j) = ^ u( 7 i/ - 7i) = Wi/ {Pd)- 

j=i j=i 

Since 71 ,, jm are conjugate and pd G X[X, X], (pd) = v^., (pd) for 1 < j < m. Recall that 
Vg{pd) = mini<j<mU 7 ,(pd). So Vg{pd) = Vg{pd), as desired. □ 

(i) 

In the general case, when the coefficients are not all different, the truncation has to be 
done iteratively. We describe a recursive process to obtain pd, the numerator of the integral 
basis of degree d in X. 

Let Pd G VxlY] be as in Lemma 7.5. Let tk G Q he the smallest exponent such that the 
truncations 

7^°^ = a-jx*! H-h aix*^, < ... < 4 

of the expansions 7 ^, 1 < j < m, are pairwise different. We truncate the expansions jj, 
1 < j < m, to degree tk-i'. 

7^ ^ = a{X*i + • • • + 

For the recursion, we define go = U]Li(Y - 7 j) and 5-0 = ( X — 7 ^^^^ • Since 4 was 

the smallest integer for which all the truncated expansions were different, the expansions 
1 < j < m, can now be grouped into sets of identical expansions, each set having the same 
number of elements. Denote by 771 ,... , 7 ^ the mutually distinct expansions, and set gi = {Y — 
r/i) • • • (X — gr) G K[X, X]. By construction go = with ui = m/r G N. 

For simplicity, we explain first how to compute recursively the element of degree m — 1, 
assuming pm-i = (X - 72 ) • • • (X - 7 ^). 

We start the Tth step by applying the whole procedure inductively to gi-i, computing 
gi and Ui such that gi-i = gfi and gWI comes from truncating the expansions of gi-i. In each 
step the degree r* of gi is smaller or equal to the degree rj_i of gi-i, and it will be equal to 1 
after a hnite number w of steps (bounded by the degree t^ of the expansions in go). For that 
value re, = I and all the expansions in g^^ are equal. The desired polynomial is 

= g^-^g^-^ • • • gfffi-^ G X[X, X]. 

We thus obtain Algorithm 6 . 

Lemma 7.6. With notation as above, let 

Pm-i = TruncatedFactor({ 7 i,.. .,7m})- 

Then Pm-i has maximal valuation at g over all monie polynomials of degree m — 1 inY. 
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Algorithm 6 Truncated Factor 

Input: A = { 7 j = + • • • + a conjugacy class of Puiseux series of finite 

length. 

Output: q G K[X, Y] of degree m — 1 in T such that v-y-^ (q) = (q), with q = (T — 72 ) • • • (T — 

7m)- 

1: Let r]i,... ,r]r he the different expansions in the set • ■ • )7m*'‘“^}) the truncations up 

to degree 

2 : p= {Y -r]i)---{Y - rir) 

3: u = mir 

4: if r > 1 then 

5: p' = TruncatedFactor({r/i,..., 77 ,,}) 

6 : return q=p^~^p'. 

7: else 

8 : return q=p^~^. 


Proof. As in the proof of Lemma 7.5, it is enough to show that u-yj(pm-i) = 1171 (Pm-i)- Let 
Y 2 -, ■ ■ ■ ,7m be the Puiseux expansions of pm-i, corresponding to truncations of the expansions 
72, • • ■, 7m of g. By construction, 11(71 - ^) = 11(71 - 7*) for i = 2,..., m. Hence Vg{pm-i) = 
V'yiiPm-i) = t’7i(Pm-i) = Vg{pm-i) as Wanted. □ 

Example 7.7. Returning to Example 4.2, the singular parts of the Puiseux expansions are 

7)°) = 7X^/2 + (-l/2i - 1/2)X^/'^ + l/4iX2 

7^ = 7X^/2 + (-l/2i - 1/2)X^/^ - l/4iX2 

7^ = 7X^/2 + (l/2i + 1/2)X^/^ + 1/47X2 

7^^ = 7X^/2 (1/27 + 1/2)X^/^ - 1/47X2 

7^°) = -7X^/2 + (1/27 - 1/2)X^/^ + 1/47X2 

7^ = -7X^/2 + (1/27 - 1/2)X^/^ - 1/47X2 

7^ = -7X^/2 + (-1/27 + 1/2)X’'/^ + l/4iX2 

7^°) = -ix3/2 + (-1/27 + 1/2)X^/^ - l/4iX2 

with 72 = —1. 

Truncating 7^°^ to degree 7/4 we obtain 

7jb = ^(1) = ix^/2 ^ (_i/2i - 1/2)X^/^ 

= 7?^ = + (1/2* + 1/2)X^/^ 

7^'^ = 7^'^ = + (1/2* - 1/2)X7/" 

7{b = = _ix3/2 (_i/2i + l/2)X^/‘^ 

hence ui = 2 and 

gi = {Y- j[\y - 7('))(y - 7W)(y - 7«) 

= + 2X^y2 + 2X^y + x® + i/4x^. 


Applying the whole procedure inductively to gi we obtain g 2 = Y'^ + X^, U 2 = 2 and gs = Y, 
773 = 2. Combining the factors, we get 

g = + 2X=^y2 + 2X^y + X^ + ^xA {Y^ + X^)Y. 
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For computing the elements of any degree d, l<d<m — 1, we can easily extend the above 
construction, leading to Algorithm 7. 


Algorithm 7 Truncated Factor General 

Input: A = { 7 * = + • • • + a '^'^a conjugacy class of Puiseux series of finite 

length; d G N, d < m. 

Output: p G K[X,Y] of T-degree d such that Vf^{p) = Vf^{p), with p the element in VxiY] 
of degree d with maximal valuation at /a- 
1 : Set rji,... ,r]r the different expansions in the set 
2 : u = [d/r\, d' = d — ur 
3: if d' > 0 then 

4 : gi = TruncatedFactorGeneral({r/i,. .., r]r}, d') 

5: else 

6 : gi = l 

7: if u > 0 then 

8 : 9= (y-Vr) 

9: return p = g^gi. 

10: else 

11: return p = gi- 


Lemma 7.8. Let g G iF[[A]][y] he a Weierstrass polynomial of degree m in Y. Let F = 
{ 71 ,..., 7 m} he the expansions of g at the origin, which correspond all to the same conjugacy 
class. Then, for any d, 1 < d < m, the output pd = TruncatedFactorGeneral(F, d) is a 
polynomial with maximal valuation at the origin in the ring A[[A]][y]/(g') among all polynomials 
in K[X, Y] monic of degree d in Y. 

Proof. Let 7 !*^^ ... , 7 m^ be the truncations of the Puiseux expansions of g up to order tk and 
9 o^ 90 : ■ ■ ■ 1 9 w— 1 t 9w—^i 9w as defined before. 

We have noted in Section 4.8 that a polynomial pd satisfying the requirements of the lemma 
can be chosen so that all the Puiseux expansions of pd at the origin are truncations of the 
expansions in F. This implies that we can take pd to be a product pd = gi^ . • • 5 ^,“ of the 
polynomials gi with appropriate exponents. To find the exponents, we note that for all i the 
polynomials giffi and gi have the same degree, but the valuation of gi at g is larger than the 
valuation of g^ff^ (since the expansions are developed up to a larger degree). Hence, to construct 
Pd, we must first take di as large as possible. Then maximize d 2 and so on iteratively. This is 
done by Algorithm 7. □ 

Setting Pd = TruncatedFactorGeneral(T, d), we compute o{F, d) by the formula 

o{F,d) = ^u(7-7), 

r]£N 

where N = { 71 ,... are the expansions appearing in pd and 7 G T. (For any expansion 7 G T 
the result of the sum is the same, because conjugating the above expression does not modify 
N.) 

Example 7.9. We carry on Example 7.7, computing all the numerators of the elements of the 
integral basis. We have obtained that the element of the integral basis of degree m — 1 = 7 is 
the product py = gig 2 g 3 , where gi, g 2 and gs have degrees 4, 2 and 1 respectively. To obtain 
the numerators of the elements of the integral basis of smaller degree d, 1 < d < 6, following 
Algorithm 7, we have to hrst take the largest possible power of gi so that the total degree is 
smaller than or equal to d, then choose the power of 52 in the same way and finally the power of 
53 . We get the following elements pe = gig2, P5 = 9193, Pi = 9i, P3 = 9293, P2 = 92 and pi = 53 . 
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The denominators are powers of x. To obtain the exponents, we compute o{r, d) for 1 < 
d < 7 hy looking at the expansions corresponding to each gi, i = 1,2,3, given in Example 
7.7. Setting Ng^ the expansions appearing in g^, i = 1,2,3, we have Yhri&N ^(7 ~ v) = 27/4, 
^^(7 -V) = 13/4 and ^^(7 - v) = 3/2 for any 7 e T. Hence o(r, 1) = 3/2, 

o(r,2) = 13/4, o(r,3) = 13/4 + 3/2 = 19/4, o(r,4) = 27/4, o(r,5) = 27/4 + 3/2 = 33/4, 
olr,6) = 271A + 13/4 = 10 and o(T, 7) = 27/4 + 13/4 + 3/2 = 23/2. The exponents in the 
denominators are the integer part of these valuations, hence the integral basis is 

f 53 92 9293 9i 9193 9192 919293 \ 
lx’x3’ x4 ’x6’ x8 ’ xio ’ xll / ’ 

7.6. Merging the integral bases for the branches. We have shown how to compute an 
integral basis when g G iir[[X]][y] is an irreducible Weierstrass polynomial. For the general case 
when g is not irreducible, theoretically we can combine all the integral bases of the branches 
following Corollary 5.10. However, as we discussed in Section 7.1, the use of the extended 
CCD is not practical. The following result, which extends Corollary 5.10, provides a different 
strategy replacing the Bezout coefficients from the Euclidean algorithm by more simple and easy 
to calculate coefficients in K[X,Y]. These coefficients as well as their integrality exponents can 
be computed based only on the singular part of the Puiseux expansions of /. 


Proposition 7.10. Let A = K[C] = K[X,Y]/(f). Let f = /o/i • • •/r be the factorization 
of f in 7r[[X]][y], where fo G i^[[y]][y] is a unit in y]], and fi,...,fr are irreducible 

Weierstrass polynomials in 77[[y]][y]. For each i = 1,... ,r, let 




(0 

Pnii-l 
Ji) 


be an integral basis as in Lemma 5.3 for the normalization of iP[[y]][y]/(/i). Let hi = 
fj. Assume fdi G K[X,Y], 1 <i <r, are such that the order at the origin of Pthi in 
-^[[^']][E/(/*) is an integer a > 0. Let B = K[[X]][Y]/{fi ■ ■ ■ fr) and set 


jSii) = 


fdihi Idihipl 


(^) 




^|cB. 

X 




Then U ... U B^^'> is a K[[X]]-module basis of B. 

Proof. For fixed 1 < * < r, by the Bezout identity on fi and fdihi there exist 9,12 G X[[X]][y] and 
e G Z>o such that gfi + n/dihi = X®. Setting ei = ord(zz) in X[[X]][y]/(/j), we have e = ei + Cj 
and ei is integer by the hypothesis on /3j. Hence v jX^^ and fdihijX^' are integral elements over 
K[[X]][Y]/{n). Call gi = v/X^T Then = big in X((X))[y]/(/,), for 1 < j < r. 

Hence, as in Proposition 5.9, the well dehned map of X[[X]]-modules 


{ti mod fi,... ,tr mod fr) 1 — E 9i 

i=l 


(dih 


-ti mod /i... fr, 


maps 0[=i X[[X]][y]/(/j) isomorphically to X[[X]][y]/(/i ■ ... ■ fr 
For 1 < z < r, let 


nil) I Pib-i jdihiPg 

W-{ 9i^pz7,9i 


(7 




, (0 
XG+ei 


I 9i 


^ihiP^Xi 


X 




We first show that as X[[X]] -modules over X[[X]][y]/(/j), + (3^)^ Clearly + 

(£(*)), because gi^^ = 1. To see that the ^ take a G Then, since gi 

is integral, giO G and a = ^^gia G proving the isomorphism. It is trivial that 

^C(®)^ + (HW) as X[[X]] modules over K[[X]]\Y]/ {fj), j / i. Hence (C^®)) + as X[[X]] 

modules over X[[X]][y]/(/i • • • fr). Therefore B^^l U • • • U is also a X[[X]]-module basis of 
B. □ 
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We address now the computation of coefficients /3i, 1 < i < t, with the required property that 
the order of Pihi at ii'[[X]][y]/ (/j) is integer. If ord(/ij) is integer, we take /3j = 1. In the general 
case, to find f3i we could use a power of hi, since ord(/i™) = mord(/ij), and we can choose m so 
that mord(/ij) is integer. (In fact, we would replace hi by hi = n7eA(^ “ 7) where A are the 
Puiseux expansions of hi and for 7 G A, 7 is the singular part of 7, so that we get a polynomial 
in K[X,Y] with the same order as hi at [y]/(/j).) Usually, however, it is more efficient 

to choose /3j as a product of the factors given by Algorithm 7. We can proceed algorithmically 
as in Algorithm 8. 

Algorithm 8 Merge Coefficients 

Input: Ai, ..., Ar, the sets of singular parts of the Puiseux expansions of the conjugacy classes 
corresponding to the factors fi,..., fr of /. 

Output: {(/Si, Ci)}i<j<,., Pi G K[X,Y] and Ci G Z>o, such that the order of Pihi in 
^[[^]W]/ifi) at the origin is Ci, where hi = Oj^i fj- 
1: for i = 1,... ,r do 
2: if ord(/ii) G Z then 

3: ,Si = 1, Ci = ord(/ii). 

4: else 

5: For each 1 < j < r, j 7^ i, and each 1 < k < dj {dj the degree of fj) set 

fj^k = TruncatedFactorGeneral(Aj, k) 

if k < dj and fj^dj = 075A “ 7 )> where 7 is the singular part of 7 as before. 

6: For each prime divisor a of the denominator of ord(/ii) G Q, take p a polynomial of 

smallest degree in Y among all the computed polynomials such that the denominator of 
ord(p) is a multiple of a (note that such polynomials always exist since hi is a product 
of some of these polynomials). 

7: Take Pi the product of these factors to appropriate powers. The exponents can be 

found by solving a linear congruence equation, choosing the solution that minimizes the 
T-degree of Pi. 

8: Ci = OT(i{Pihi). 

9: return {{Pi,Ci)}i<i<r. 


To merge the integral basis from the branches applying Proposition 7.10, it remains to truncate 
the elements hi,\ < i < r, to polynomials in K[X,Y]. Note that the coefficients Cj, 1 < z < r, 
can be computed from the singular part of the Puiseux expansions of /. If Cc is the maximum 
order of the coefficients Ci, then we know that for any polynomial appearing in the construction 
of the integral basis the integrality exponent will be at most Cc + E{f). Hence, by Remark 5.4, 
we can truncate all the numerators to degree Bc + E{f) in X. We obtain Algorithm 9. 

Remark 7.11. To speed up the computation of the integral basis, we first compute the order e 
of y in K[X, Y]/{g) and add the elements , 0 < z < m = degy(y), to since 

those simple elements also belong to the normalization. This is an improvement over Remark 

5.4. 

Example 7.12. Let A = KIX,Y]/ ((Y^ + X^)(Y^ - X^) + Y^) = K[x,y] as in Examples 4.5 

and 5.11. There are two conjugacy classes of expansions at the origin. Pi = {71,72,73} and 

r2 = {74,75}. We apply Algorithm 8 to compute Pi and hence as in Proposition 7.10. 

We have fi = (Y — 7i)(y — 72)(H — 73) and hi = {Y — 74)(y — 75). The integral basis of 

2 

iL[[A]][y]/(/i) is {l,y, ^}. Evaluating the expansions from Ti in hi, we see that the order of 
hi at fi is 4/3. Applying Algorithm 7 to hi we get TruncatedFactorGeneral(/ii, 1) = Y. The 
order of y at fi is 2/3. Hence yhi has order 2 at /i, which is integer. So we can use Pi = Y. 
We get 

„(i) _ / yhi yhiy yhiy‘^\ 

\ ’ x3 / ■ 
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Algorithm 9 Integral basis for Weierstrass polynomial 


Input: L = {{A, /i},..., {A, fr}}, where A = , 7i,mJ is the set of singular parts of 

the i-th conjugacy class of expansions that vanish at the origin of a polynomial / G K[X, y] 
monic in Y and fi is the corresponding factor of /; {(/3i, ci(A; Cr)}, the output of 
MergeCoef f icients({A, • • ■) A})- We assume fi, 1 < i < developed up to X-degree 
ec + E{f), with Cc = maxi<j<r c*. 

Output: {{po,eo ),..., {pm-iXm-i)}, Pi £ -^[W][y] monic in Y of degree i and e* G Z<o, such 
that • • •, A } is an integral basis for the normalization of A[[X]][y]/(/i • • • A). 

1 : m = degy(/i---/r) 

2: for 7 = 1,..., r do 

3: ^i = Y\j^ifj 

4: for d = 0,..., m* — 1 do 

5: qd = TruncatedFactorGeneral(A, d) 

6: = the integrality exponent of qd in K[X,Y]/{fi). 

7: pd = bihiqd, ed = Ci + e{qd) 

o. r7{i) _ f _P0_ _P1_ Prru-l \ 

o- ’ ■ ■ ■ ’ / 

9: From compute the integral basis {po,Pi/x^^, ■ 

in the proof of Lemma 5.3. 

10: return {(po, eo), ■ • •, (Pm-i, em-i)}. 


,Pm-i/x^^-^,}, as indicated 


This is more simple than using 6i = — 2X^y^ — 3X^ — 2X^y + XY^ — Y^ — Y as in Example 

5.11. 

We next consider a more complicated example where the coefficient A has larger degree. 

Example 7.13. Let f{X,Y) = {Y^ - 6XW^ - 2X’^Y‘^ + UX^Y"^ - UX^^Y - 8X^){Y^ - 
2YX^ — 2X^)(y^ + X’^) + X^° and A = K[X,Y ]/(/). The Puiseux expansions of / are 

71 = riX3/2 + x7/3 ^ ^ ^ ^ ^^^7/3 ^ ^ ^ 

72 = riX3/2 + rgX^/^ ^ ^ ^ ^^^3/2 + 7^3 ^ ^ 

73 = riX3/2 + r4X^/3 ^ ^ ^ ^^^3/2 + 7^3 ^ ^ 

74 = r2X3/2 + x7/3 + , . . ^ = ^^7^7/2 ^ ^ ^ 

75 = r2X3/2 + rsX^/^ ^ ^ ^ ^^7^7/2 ^ ^ _ 

with ri, r 2 the roots of — 2 = 0, ra, r 4 the roots of + A +1 = 0 and , rg the roots of A^ +1 = 
0. They correspond to three conjugacy classes Ai = { 71 , 7273 , 74 , 75 , 76 }; ^2 = {77,78} and 
^3 = { 79 , 710 }- We show how to compute A- Let /i, A and /a be the Weierstrass polynomials 
corresponding to each conjugacy class and hi = A/a- 

The order of hi at A is 41/6 (recall that the order of hi at A is the order in X of (7 — 
77 )... (7 — 710 ) for any expansion 7 of hi). Applying Algorithm 7 to A 2 and A 3 for c = 
1,2 we get the polynomials {y,y^ — 2yx^ — 2x^'\ and {y,y‘^ + x'^}. We have ord/Ad) = 3/2, 
ord/j(d^ “ 2yx^ — 2x^) = 23/6, ord/Ad^ + x'^) = 3. 

We can take /3i = y®Ad^ “ 2yx^ — 2x^)®2 for appropriate ei,e 2 G Nq. The exponents must 
satisfy the equation ei| + 62 ^ + ^ G Z. The corresponding linear congruence equation is 
9ei + 2362 + 41 = 0(6). The solution that minimizes the T-degree of A is ei = 1 and 62 = 2. 
We get A = d(d^ “ 2yx^ — 2x^)^. 


7.7. Ad-hoc Algorithm for the Case of oue Characteristic Expoueut. We describe an 
algorithm, which allows us to write down an integral basis directly in the case of a singularity 
with only one conjugacy class of Puiseux expansions. 
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Lemma 7.14. Let f G -f^[[2;]][y] be an irreducible Weierstrass polynomial with respect to y and 
degy f = n. Let y{x) be a Puiseux expansion, and y{x) = d, m > n and 

gcd(m,n) < n. Let ko = n, ki = m, k 2 ,---,kg be the characteristic exponents and let e be 
a primitive n-th root of unity. The following holds (cf. [9, Lemma 5.2.18(1)] and the proof 
thereof): 

n 

(1) f = Yliv-Vi^'x)) 

i=l 

(2) For j = 1,... , g denote by Nj the set of all i G {1, ... ,n} such that 

h I • , ko 


gcd(A:o,.. ■,kj-i) 


i and 


Then 

for all i G Nj. In particular, if g = 1 then 


gcd(A:o, ...,kj 


kn 


\i. 


oi:dj:{y{x) - y(e*x)) = N. 

n 


ki 

oidx{y{x) - y(e*x)) = — 

n 


if i is not a multiple of n. 
(3) We have 


d f k ■ 

orda: -^{x, y{x)) = } (gcd(A;o,..., kj-i) - gcd(/co, • ■ •, kj)) 
ay ^ n 

j=i 


= Int 




for all i. 


Proposition 7.15. With notation as above we have: 

d£ 

(1) For e = orda; ^{x,y{x)) the element is integral over i^[[x]] and e is maximal. 

(2) Let 


e = 


dn-lf 


1 — 1 


Then kLN — jg integral over /C[[x]], e = 
(3) If g = 1 then 


and e is maximal. 


d"-^f ^ 

dy'^~^ dy 

, . . . , 


X^ 


X^ 


with 


e,; = 


dn-if 


form an integral basis of K[[x,y]]/ (/) over K[[x]]. 

We now prove Proposition 7.15. 

Proof. Choose LI C {1,..., n} with jil] = d and Into maximal. Then 

P-=Y\{y- 

j&n 

is a polynomial of degree d with respect to y and OTdxp{x,y{x)) is maximal. Let 

e = [ovdxp{x,y{x))\ . 

By Lemma 7.8, for some approximation p of p, we can choose degree d element 

in the integral basis. 
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We obtain (1) for d = n — 1 and (2) for d = 1, and Into is independent of fi. The same holds 
true for (3) in case g = 1. □ 

Remark 7.16. If/ = y^ — 2x‘^y^— ^x^^y + x^— then y(x) =xi+x^ is a Puiseux expansion, 

g = 2 and (3) in Proposition 7.15 does not hold. 


Proof. We compute 


orda; 


dfl 

dy^ 


{x,y{x)) 


= 4. However, 


p={y- y{-x)) {y - y{ix)) 


gives [orda;p(x, y(x))J =6. 


□ 


7.8. Computation of the Local Contribution to the Integral Basis. Finally, assuming 
that P = {X, Y) is the only singularity at X = 0, we compute the local contribution to the 
integral basis from the integral basis of iF[[W]] [y]/(/i ■ ■ ■ fr) using Proposition 6.1. We describe 
the complete procedure in Algorithm 10. 


Algorithm 10 Local contribution to the integral basis 

Input: / G K[X,Y\ irreducible polynomial, monic of T-degree n, with only one singularity at 
X = 0, located at the origin. 

Output: A basis of the minimal local contribution of K[X,Y]/(f) at the origin. 

1: Compute Ao the set of the singular parts of the Puiseux expansions of / at A = 0 that do 
not vanish at the origin and Ai,..., A,, the sets of singular parts of the conjugacy classes of 
Puiseux expansions of / at A = 0 that vanish at the origin. 

2: Compute E{f) as indicated in Section 4.8. 

3: Compute {(/3i, Ci)}i<i<r = MergeCoeff icients({Ai,..., A^}). 

4: €c — maX]^<^<y. Ci- 

5: {/o, = Splitting(/, E{f) + Cc), where /o corresponds to Aq and /i,..., /^ cor¬ 

respond to Ai,..., Aj,. 

6: mo = deg(/o), m = n - rriQ 

7: ..., = IntegralBasisForWeierstrassPolynomial({(Ai,/i), ..., 

(A^, fr)},{i/3l, Cl), . . . , (Pr, Cr)}) 

8: for z = 0 ,..., deg(/o) - 1 do 

9: Pi = y\ Cj = 1 

10: for z = 0 ,..., m — 1 do 
11- PrriQ+i — fo ' Pi^ Cmo+i = 6 (*) 

12: return B = {pq/x^^^,... ,pn-i/x^'^-^}. 


To compute a (global) integral basis of A over K[x] we can now use Proposition 3.1. 

Remark 7.17. In the presence of conjugated singularities, to get a better performance, our 
local algorithm can handle groups of conjugate singularities simultaneously, in a similar way as 
in [21, Section 4]. If I C A [A, A] is an associated prime of the singular locus, corresponding to a 
group of conjugate singularities, we apply a linear coordinate change if necessary, so that no two 
of these singularities have the same A-coordinate. Then we can find polynomials qi,q 2 £ A[A] 
such that I = {qi{X),Y — q2{X)). We take a a root of qi{X) and translate the singularity 
(a, 52 ( 0 :)) to the origin. We compute the local contribution to integral basis at the origin and 
apply the inverse translation to the output. The common denominator of the resulting generators 
will be a power of x — a. We replace {x — a) by qi{x) in the denominators and we eliminate 
a from the numerators by considering a as a new variable and reducing each numerator by 
the numerators of smaller degree (written all with the same common denominator), using an 
elimination ordering a ^ y ^ x. Since an integral basis over the original ring always exists, the 
elimination process is guaranteed to eliminiate a from the numerators. 
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Example 7.18. Let A = K[X,Y] and f{X,Y) = Y^ — {X‘^ — 2)^. The singular locus contains 
only one primary component {X"^ — 2,Y‘^), with radical {X‘^ — 2,Y). It consists of the two 
conjugated points (—\/2,0) and (\/2, 0). We take a = \f2 and compute the local contribution 
at (a, 0) translating that point to the origin. After the inverse translation, we get the integral 

basis of the local contribution |l,y, 

The local contribution to the integral basis at the conjugate singularity is |l, y, Hence 

the global integral basis is |l,y, simple case, we did not need to eliminate a 

from the numerator.) 

Example 7.19. Let A = K[X,Y] and f{X,Y) = (Y - Xf - (X^ - 2)^. Now the radical of 
the singular locus is the prime ideal (X^ — 2, Y — X). It consists of the two conjugated points 
{—y/2,—y/2) and {V2,\/2). We take a = y/2 and compute the local contribution at {a, a). We 
get the integral basis of the local contribution 

-2ay +2 } 

X — a J 

To eliminate a from the last numerator, we write all the fractions with the same denominator 
f y{x-a) y - 2 ay+ 2 \ j ^ reduce the last one to get |l,?/, ^ -2x^+2 ! Hence 

I X—a ’ x—cx. ’ x—a J’ G I x—cx. J 

the global integral basis is |l,y, ^ 



8. Timings 

We present timings, comparing the implementation of our integral basis algorithm^ in SINGU¬ 
LAR with obtaining an integral basis via the local normalization algorithm^ outlined in Section 3, 
with the implementation of van Hoeij’s algorithm'^ in Maple [18] and with the implementation 
of the variant of the Round 2 algorithm® in Magma [6, 14]. 

We compute integral bases for A = Q[Ai, Y]/ {f) with polynomials / as specified. All timings 
are in seconds, taken on an AMD Opteron 6174 machine with 48 cores, 2.2GHz, and 128GB of 
RAM running a Linux operating system. A dash indicates that the computation did not hnish 
within 6000 seconds. We only use parallel computations for the decomposition of the singular 
locus. The parallelization of the integral basis algorithm and a modular approach following 
the strategy of [3] is subject to ongoing work. Recall that for obtaining the integral bases, 
singularities at inhnity of the curve {/ = 0} do not matter. 

8.1. One Singularity of Type A. The plane curves with defining equation f{X,Y) = Y^ + 
jjAfc+i _|_ fc > 1, d > 3 have exactly one singularity at the origin, which is of type A^. 


k 

d 

SiNC 

INTBAS 

rULAR 

NORMAL 

Maple 

Magma 

5 

10 

0 

0 

0 

0 

5 

100 

0 

0 

1 

168 

5 

500 

0 

1 

49 

* 

50 

60 

0 

0 

1 

294 

50 

100 

0 

1 

2 

10751 

50 

500 

0 

0 

76 

* 

90 

100 

0 

1 

4 

* 

90 

500 

0 

1 

102 

* 

400 

500 

0 

3 

346 

* 


^column Singular intbas in the tables 
^column Singular normal in the tables 
■^column Maple in the tables 
^column Magma in the tables 
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8.2. One Singularity of Type D. The plane curves with defining equation/(X, T) = 
k > 3, d > 3 have exactly one D^+i-singularity at the origin. 


k 

d 

SiNC 

INTBAS 

,ULAR 

NORMAL 

Maple 

Magma 

5 

10 

0 

0 

0 

0 

5 

100 

0 

1 

1 

1683 

5 

500 

3 

0 

53 

* 

50 

60 

0 

1 

6 

312 

50 

100 

0 

1 

17 

3480 

50 

500 

3 

1 

808 

* 

90 

100 

0 

1 

51 

* 

90 

500 

3 

1 

103 

* 

400 

500 

3 

4 

2326 

* 


8.3. Ordinary Multiple Points. We consider random curves of degree d with an ordinary 
fe-fold point at the origin. The defining polynomials were generated by the function polyDK 
from the SINGULAR library integralbasis. lib (using the random seed 1231). 


k 

d 

Singular 

Maple 

Magma 



INTBAS 

NORMAL 



5 

10 

0 

2 

0 

0 

15 

20 

0 

7784 

1 

4 

15 

30 

1 

* 

28 

124 

20 

25 

0 

* 

2 

18 

20 

30 

1 

* 

19 

42 


8.4. Curves With Many Singularities of Type A. The plane curves with defining equations 

/5,, = + 2(x”y” - x”y” + y’"z") 

were given in [7] and have 3n singularities of type A^-i if n is odd. We substitute Z = X —2y+1. 


n 

Singular 

Maple 

Magma 


INTBAS 

NORMAL 



5 

0 

1249 

1 

1 

7 

1 

* 

7 

8 

9 

18 

* 

30 

59 

11 

56 

* 

231 

251 


8.5. More General Singularities. We now consider some examples of curves which have 
singularities of a type other than ADE or ordinary multiple points: 

(1) / = + 2lX^* - SX^^Y + 6X^3 + 16Xi2y _ 20X^^Y^ + X^^ - SX^V + 36X^°Y^ - 24X®W - AX^Y^ + 

16X®Y® — 26X^Y* + 6X®y^ — 8X®y® — 4X^y® + y®: one singularity at the origin with multiplicity m = 8 and delta 
invariant 6 = 42, a node, and a set of 6 conjugate nodes. 

(2) / = (y*+ 2 X®y 2 + x®+ x5y)® yxiiyiU one singularity at the origin with m = \2 and 5 = 133. 

(3) / = (y5 + Yix’^ + 2X®)(y3 + 7X*)(Y’^ + 2Xl2)(yll + 2X^^) + y^O; one singularity at the origin with m = 26 and 
S = 523. 

(4) / = (yl5 + 2X3®)(yl9 + 7X52) +y36. one singularity at the origin with m = 34 and (5 = 1440. 

(5) / = (yl5 _|_ 2X5®)(yl2 + 7X52) _|_ yioo. higher degree, but same type of singularity. 

(6) / = y^® + Xy^5 + X^y5 + X5 + 2X'* + X®: one double point with 5 = 2 and one triple point with 5 = 19 (see [21, 
Section 6.1]). 

(7) / = y200 _|_ XY^^ + X^y5 + x5 + 2X^ + X®: higher degree, but same type of singularity. 

(8) / = (y35 + y34x7 + 2X38)(y33 + 7X44)(y37 + 2X52) yllO. one singularity at the origin with m = 105 and 
S = 6528. 
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No. 

y-degree 

SiNC 

INTBAS 

,ULAR 

NORMAL 

Maple 

Magma 

(1) 

8 

0 

* 

0 

0 

(2) 

12 

9 

* 

1 

1 

(3) 

30 

11 

* 

5 

31 

(4) 

36 

1 

* 

5 

59 

(5) 

100 

1 

* 

34 

* 

(6) 

40 

0 

2 

1 

9 

(7) 

200 

1 

3 

12 


(8) 

no 

* 

* 

* 



In Example (8), Singular and Maple do not finish due to the computation of the decom¬ 
position of the singular locus of the curve. See Section 8.6 for the timings of the computation 
of the local contribution to the integral basis at the origin. 

We note that in most cases, our proposed algorithm is much faster than the algorithms 
implemented in Maple and Magma. 

8.6. Detailed Analysis of Some Examples. The computation of an integral basis with our 
algorithm has two major components. First, we decompose the singular locus into associated 
primes, secondly we compute the local contribution to the integral basis at each associated prime 
and then we merge the results. 

In Maple a similar strategy is followed, computing first all the x-coordinates of the singular 
points. This first step is part of both integral basis approaches, and can be time in consuming in 
some examples. To analyze in more detail the difference between the integral basis computations, 
we provide here timings for the computation of the integral basis at the origin for some examples 
where the origin is the only singularity. This can be specified in Singular and Maple by certain 
input options. 


Example 

k 

d 

Singular 

INTBAS 

Maple 

8.1 

5 

500 

0 

49 

8.1 

50 

500 

0 

75 

8.1 

400 

500 

0 

346 

8.2 

5 

500 

2 

49 

8.2 

50 

500 

2 

571 

8.2 

400 

500 

2 

1575 

8.3 

15 

30 

0 

26 

8.3 

20 

25 

0 

2 

8.3 

20 

30 

1 

17 

8.5 (2) 



9 

1 

8.5 (3) 



4 

2 

8.5 (5) 



1 

16 

8.5 (8) 



25 

1483 


We observe that for the examples in Section 8.1 the time required for the singular locus 
decomposition is not significant. For the examples in Section 8.2, Maple uses a significant 
amount of time for this task, but still the most time consuming part is the computation of 
the integral basis at the origin. For the examples in Section 8.3, in our implementation most 
of the time is used for computing the decomposition of the singular locus. The computation 
of the integral basis at the origin is significantly faster than in Maple. For the examples in 
Section 8.5 various situations occur. In Example 8.5 (2) the time for the initial decomposition 
was not signihcant, and computation of the integral basis at the origin in our implementation 















INTEGRAL BASES VIA LOCALIZATION 


33 


is slower than Maple. In this example, the algorithm runs into an algebraic field extension of 
high degree. At current state, the handling of such extensions in Singular is not optimal. 
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